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THE INITIAL IONIZATION PRODUCED FISSION FRAGMENTS! 


Abstract 


The ionization produced individual fission fragments over various lengths 
track starting from the origin was measured pulse ionization chamber 
provided with grid. The distribution electron pulse sizes for the complete 
track showed the well known two peaks. the length track was diminished, 
the peaks corresponding the high and low energy groups fragments 
approached each other and merged together track length about mm. 
nitrogen. smaller distances measured from the beginning the track the 
difference ionization the two groups fragments was found small, 
and the two groups could not resolved. From the progressive change shape 
the pulse spectrum tentatively concluded that the low energy group 
produces the greater initial ionization. The average energy spent ionization 
the fragments the first 1.1, 2.3, and 5.3 mm. their tracks was measured 
nitrogen and found 7.6, 15.0, and 34.6 Mev. course 
this work effect due the formation negative ions electron attach- 
ment was observed the chamber filled with argon. Examination this 
effect suggests that, high collecting fields, formation negative ions, perhaps 
impurities such oxygen and water vapor, more pronounced argon than 
nitrogen. 


Introduction 


The general features the loss energy fission fragments passing 
through matter have been described Bohr (2), who emphasized the essential 
differences mechanisms which and fission fragments 
lose energy. the first place, the effective charge essen- 
tially equal the nuclear charge during most its range and decreases 
only towards the end the range. the contrary, the effective charge 
fission fragment much less than the nuclear charge and decreases continu- 
ously from the beginning the track. Consequently, the specific ionization 
fission fragment can compared that between the 
maximum the Bragg curve and the end its range. energy region 
which the specific ionization increases with decreasing velocity for particles 
massive fission fragments should occur energies much higher than 
those actually attained fission, while the case this region 
occurs energies above one Mev. Thus, for the maximum value 
the specific ionization occurs not far from the end the range, while for 

fission fragments actually occurs the beginning the tracks. 


Manuscript received January 1948. 
Contribution from the Nuclear Physics Branch, Chalk River Laboratory, Atomic Energy 
Division the National Research Council Canada. Issued N.R.C. No. 1745. This paper 
was presented September 1946, the Nuclear Physics Conference Montreal, Que. 
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There another mechanism which charged particles lose energy 
passing through matter: the ‘indirect’ ionization nuclear collisions. This 
effect much greater importance for fission fragments than for 
fact, energy loss nuclear collisions, which occurs with significant 
frequency only the very end the track takes place all 
along the track fission fragment. However, even for fission fragment, 
the amount energy lost this process small comparison the energy 
lost direct ionization, except the very end the track. 

Calculations the energy loss fission fragments passing through 
matter have been made Lamb (11), Bohr (2), and Knipp and Teller (10). 
The detailed comparison these calculations with experiment is, however, 
made difficult the wide distribution mass and energy the fragments. 
The energy loss has been investigated cloud chamber 
and Lauritsen (1). They used the recoil tracks that occur along the track 
the fission fragment deduce its velocity various points. Their method 
most accurate towards the end the range where recoils are most frequent. 
Information about the initial part the range was obtained Lassen, who 
performed magnetic deflection experiments and deduced the initial effective 
charges the light and heavy groups fragments (12) and the effective 
charge and energy the light group fragments after they had traversed 
various thicknesses mica (13). Lassen obtained results with which Bohr’s 
theory does not agree detail. 

view this disagreement and the lack more complete experimental 
data about the initial part the range, measurement the initial ionization 
fission fragments was undertaken. 


Experimental Method 


The principle the method follows. thin sample fissile material, 
this case uranium 233 (16), was placed the center hemisphere. The 
sample holder and the hemisphere constituted two electrodes pulse 
ionization chamber. varying the pressure the gas the chamber the 
ionization produced various distances beginning the origin the track 
fission fragment could investigated. The hemispherical arrangement 
ensures that the length track the gas the same for all 
was chosen preference parallel-plate geometry and well collimated 
beam fission fragments avoid loss ionization the channels the 
collimator. 

When the simple hemispherical chamber with two electrodes was filled 
with argon and tested, was found unsatisfactory owing the long 
collection time positive ions the electric fields that could applied. 
Mr. Wilkinson suggested form grid chamber which this difficulty 
avoided. was adopted for the experiment and shown Fig. 
hemispherical geometry retained but the support for the source made 
very small. The potential the support adjusted that electrons produced 
the ionizing fragments are not collected but move past the collecting 
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electrode. grid placed front the collecting electrode shields from the 
effect the motion the positive ions. This grid consisted grating 
parallel wires, 0.13 mm. diameter, spaced intervals 1.4 mm., and 


Fic. pulse ionization chamber. 


was placed 1.1 cm. front the collecting electrode. the 
theory grid chambers that has been given Bunemann, Cranshaw, and 
Harvey (5), the residual effect the positive ions produces variation less 
than the pulse size this chamber. 


The chamber was irradiated with neutrons the low power, heavy water 
pile. The fission pulses were amplified linear amplifier, displayed 
the screen cathode-ray oscilloscope, and photographed mm. film. 
The deflection the screen varied linearly with the voltage applied the 
input the amplifier deflections corresponding pulses 200 
gain control that consisted bank calibrated resistors constituting 
anode load was used the was placed the first stage ensure 
linearity all gain settings. The bandwidth the amplifier extended from 
0.6 kc. per sec., which frequencies the response fell half the 
maximum. This width was sufficient ensure that the pulse height was 
independent any variation the collection time the electrons that would 
encountered the experiment. The pulses the oscilloscope screen were 
photographed moving film. The horizontal sweep was not used. The 
direct spot the stationary beam the screen was covered with narrow 
slit eliminate fogging the film along the base line. This arrangement 
gave well defined base line the photographic record from which precise 
measurements the pulse heights could made. 


Preparation the Fission Source 


Since the initial ionization the fission fragments measured, the 
loss energy the layer fissile material must kept small possible. 
idea the amount energy lost the layer was obtained calculating 


HT NEGATIVE 
SOURCE ! OLLECTING 
ay i ELECTRODE 
V7 
' 


118 CANADIAN JOURNAL RESEARCH. VOL. 26, SEC. 


the mean energy lost initially monokinetic group fission fragments 
layer thickness was assumed that the relation between specific 
energy loss and the residual range the form 


where Over the initial part the track this relation good 
approximation the results and Lauritsen The 
calculation was carried out for the case when the layer fissile material was 
covered with cylindrical collimator, since some collimation was found 
necessary. The behavior the collimator was idealized assuming 
that the number particles emerging per unit solid angle independent the 
angle emergence. Since fact the number particles per unit solid angle 
emerging from the collimator decreases almost linearly with the tangent 
the angle emergence, falling zero the maximum angle, the expression 
written for the mean energy loss will upper limit. The expression 
derived the above assumptions 


where mean energy spent the layer fragments that emerge from 


the collimator, 
initial energy the fragment, 


thickness the layer, 
range the fragment the fissile material, 
and maximum angle emergence fragments from the cylindrical 


collimator, measured from the normal the layer. 


The layer fissile material was uranium oxide weighing about 0.04 mgm. 
The collimator consisted holes, 0.75 mm. diameter, drilled 
metal sheet 0.5 mm. thick. Considering the light group fission fragments, 
which have initial energy about Mev. and range mgm. per 
uranium oxide, the mean energy loss the layer about Mev. 

The error the measured ionization the gas the chamber, caused 
the loss energy the solid layer, depends upon the pressure the gas. 
more important low pressure and sets lower limit the pressure 
which significant results can obtained. However, the effect the layer 
even low pressures not important may appear first sight. 

Consider solid layer equivalent depth gas and depth gas 
the chamber. Fission fragments that originate the bottom the layer 
lose amount energy equal that spent the first cm. the range 
before they enter the gas. Thus traversing the chamber they lose 
amount energy equal that spent between and (instead between 
Hence the energy lost the layer replaced amount equal 
the energy spent between and The net effect only small reduction 
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the energy spent the gas, this reduction being small fraction the 
energy spent the layer. The amount which the energy spent the 
gas reduced can estimated using the approximate energy-range relation 
given above. The maximum reduction energy spent the gas due 
the thickness the layer fissile material chosen was estimated about 
for chamber depth equivalent mm. air. The layer chosen (0.04 
mgm. per was therefore sufficiently thin for the measurements. 


The layer fissile material was prepared evaporating solution 
uranium nitrate pyridine platinum foil. About was 
converted nitrate and dissolved pyridine. This material was evaporated 
the deposit was ignited. The final deposit appeared have uniform 
yellow color. The gain weight the foil showed that appreciable 
impurity had accumulated from the evaporation the pyridine. small 
portion this layer was selected and used the source. 


Chamber Characteristics 


view the unusual design the chamber, detailed study was made 
the variation pulse size with the voltages applied the electrodes. 

sample polonium was placed the source electrode and the chamber 
filled with argon (99.8% pure) pressure just sufficient prevent the 
a-particles from reaching the hemisphere. The potentials the electrodes 
were adjusted give the maximum pulse size, and the distribution sizes 
was investigated. well defined peak was obtained, but addition about 
one-third the pulses were smaller sizes. These small pulses were shown 
not due the thickness the source. was thought that electrons 
released near the rim the hemisphere might collected the case the 
chamber and hence not contribute the pulse the collecting electrode. 
This explanation was shown right, since the small pulses were, fact, 
eliminated placing shallow collimator over the source. This collimator, 
described the previous section, reduced the number small pulses less 
than the limited the emission ionizing particles angle 
that indicated dotted lines Fig. 

The variation pulse size with the potentials the grid and source elec- 
trode was studied, the potential the hemisphere being kept constant. 
(The electron-collecting electrode connected ground through resistance.) 
Two the characteristics obtained are shown Fig. 


The grid characteristic Fig. has well defined plateau when its 
potential made sufficiently negative prevent its collecting any electrons. 
the upper end the plateau there falling-off mean pulse size, which 
associated with broadening the spectrum pulses. This due 
the increased ratio leakage field (between hemisphere and case) collecting 
field (between hemisphere and grid) the potential the grid made more 
negative. The effect similar that produced reducing the degree 
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collimation the source. fact, the length the plateau the grid 
characteristic was found depend the degree collimation used. one 
would expect, the plateau was shorter when wider collimator was used. 


PULSE HEIGHT 


GRID POTENTIAL 


PULSE HEIGHT 


POTENTIAL 


Fic. Graph pulse height against the grid potential when the source potential was 
0.95. 

Graph pulse height against the source potential when the grid potential was 
0.25. 


Potentials are expressed units that the hemisphere. The 
collecting electrode grounded through resistance. 


characteristic (B) for the source electrode also shown Fig. 
the potential the source made more negative, fewer electrons are collected 
and the pulse height increases steady value. The potentials source 
and grid that were chosen optimum values are indicated arrows Fig. 


After the preliminary investigation the characteristics the chamber 
with polonium was completed, the saturation voltage needed 
collect the ions produced fission fragments was investigated various 
pressures argon. The voltages were applied the electrodes means 
potentiometer that varied the total voltage without changing the fractional 
voltages the different electrodes. The pulse heights were measured 
directly graduated oscilloscope screen. The mean pulse size was deter- 
mined series applied voltages, and from the curve mean pulse size 
versus applied voltage the saturation voltage was determined. 


pressure (70 cm. mercury) sufficient stop the fragments the 
gas the chamber, saturation was obtained with —1400 applied the 
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hemisphere. interest compare this voltage with that required 
collect the ions from polonium the same pressure gas. 
that case —300 the hemisphere was required produce saturation. 


Investigation the saturation characteristics for fission fragments lower 
pressures argon revealed phenomenon that proved importance 
and led the use nitrogen instead argon. 


argon pressure cm. mercury, the mean pulse size reached 
maximum with —100 applied the hemisphere. increasing the voltage 
above 100 the mean pulse size decreased, and with the pulse size 
was only one-tenth its maximum size. similar behavior was observed 
higher pressures. check the reproducibility the phenomenon was 
made refilling the chamber cm. pressure: the effect was observed 
again but the variation the pulse size with voltage could not reproduced 
exactly. 

Observation the pulse shape showed that the decrease pulse height 
was associated with the appearance second pulse that occurred about 
one millisecond later. This delayed pulse rose much more slowly than the 
first and was consequently not amplified much. 


was apparent from these observations that some electrons were becoming 
attached gas molecules during their passage across the chamber. The 
negative ions formed move much more slowly than the electrons, and, since 
pulse produced the collecting electrode only after the ions have reached 
the grid the chamber, the slow pulse from the negative ions delayed with 
respect that due the electrons. 


The lack reproducibility the phenomenon refilling suggests that the 
negative ions are formed with impurities the argon. The argon used was 
99.8% pure; the effect, however, was also observed with argon 99.96% pure. 
The impurities were not identified, but for the purpose the present argument 
they are assumed oxygen and water vapor. For both these molecules 
the probability negative ion formation electron attachment rises rapidly 
with electron energy above certain region energy (3, 
measurements electron attachment oxygen water vapor present 
small amount argon are not available. Nevertheless, believed that 
the range electron energies attained the present chamber extended into 
the region where the probability electron attachment rising rapidly. 
This would explain the dependence the pulse size the applied electric 
field. The collecting field the region traversed fission fragments 
sufficient collect electrons, but would inadequate prevent recom- 
bination positive and negative ions. Consequently, the fact that the 
‘negative ions produced pulse rather than recombining with the positive ions 
indication that they were formed away from the positive ions, 
away from the region traversed fission fragments (see Fig. 1). 


One can thus visualize the effect follows: the velocities the electrons 
the region beyond the sample, where positive ions are not produced, increase 
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with increasing field, that the higher the field, the greater the amount 
negative ion formation from impurities. These negative ions can now 
reach the collecting electrode notwithstanding their low mobility, because 
the region considered there are positive ions produced. 

The effect described above introduces serious difficulty for the measure- 
ment the initial specific ionization fission fragments when argon used. 
The following considerations suggest that other gases may more desirable 
connection with this difficulty. 

The motion electrons gases has recently been summarized Healey 
and Reed (8). The motion under the action electric field diffusion 
process that characterized two velocities, the agitational velocity and 
the drift velocity. The drift velocity will determine the collection time 
the electrons. the mean, the direction the field, the generally 
much larger agitational velocities. will the agitational velocity that 
determines the formation negative ions, however. Now the agitational 
velocity argon much higher than most other gases for given applied 
field. Consequently the degree negative ion formation impurities 
oxygen and water vapor will considerable even for very small quantities 
impurity, as, for example, indicated the occurrence negative ion 
formation argon 99.96% purity mentioned above. 

nitrogen the agitational velocity electrons much less than argon 
for the same drift velocity. This essentially due the larger loss energy 
that occurs collision electron with nitrogen molecule. this 
account one may expect nitrogen considerably less sensitive impurity 
than argon and particular that negative ion formation will much less 
probable than argon for given drift velocity. 

The nitrogen used was the commercial grade, which was purified passing 
over heated copper remove oxygen. indication negative ion 
formation was obtained with this gas even with field strength five times 
great that which became important the case argon. The 
voltage required produce saturation the pulses was, however, 
great with argon under similar conditions. Consequently, was decided 
use nitrogen order avoid the extreme purity gas that would have 
been required had argon been selected.* 


Results 


order verify that the operation the apparatus whole (chamber 
and associated electronic equipment) was satisfactory, the distribution the 
total energy fission fragments from uranium 233 was measured pressure 
(72.5 cm. mercury) sufficient stop all the fragments the chamber. 
Unfortunately was not possible use nitrogen for this measurement 
the necessary saturation voltage could not applied the chamber owing 


Margenau al. (15) have shown that electron attachment argon takes place measurable 
degree even after extreme precautions remove have been taken. 
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electrical breakdown. With argon the plateau the characteristic 
mean pulse size versus total voltage, however, extended from 1600 2000 
(above which the decrease due negative ion formation set in). Conse- 
quently, was thought that even with argon this experiment the total 
energy fission fragments would give significant results for control purposes. 
The pulse heights were measured the photographic record the oscilloscope 
screen. The spectrum pulses obtained for the total energy fission frag- 
ments from uranium 233 shown Fig. calibration the apparatus 


ENERGY (MEV) 


NUMBER PER UNIT ENERGY INTERVAL 


ENERGY 


Fic. Energy distribution 1038 fission fragments from uranium 233 argon. 
The depth the chamber exceeded the range. 


Energy distribution 2040 fission fragments from uranium 233 nitrogen, 
made with electronic pulse The depth the chamber was equiv- 
alent 8.8 mm. path standard nitrogen. The analyzer was not calibrated 
Mev., hence the abscissae are arbitrary proportional the energy 
spent tonization. 


was obtained photographing the pulses from polonium with 
pressure argon sufficient absorb the whole track within the gas the 
chamber. The calibration was performed different setting the gain 
the amplifier. Knowing the relative values the gain the two cases the 
energy the fission pulses could determined. The energy spectrum 
the fission fragments close agreement with the energy spectra obtained 
Jentschke (9) for the fission fragments from uranium 235 and uranium 238. 
was concluded that the operation the chamber was satisfactory. 


The actual experiment the initial ionization fission fragments was 
carried out using nitrogen for the reasons already explained. separate 
calibration was performed with polonium nitrogen. order 
able apply this calibration all subsequent measurements the gain 


ro 
COCO 
q 
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the amplifier was measured each case that changes gain could 
allowed each pressure the relation between the mean pulse height and 
high voltage was investigated, and suitable collecting voltage chosen. The 
distributions obtained for the energy spent ionization the first 5.3, 2.3, 
and 1.1 mm. the tracks nitrogen reduced normal temperature and 


NUMBER PER UNIT ENERGY INTERVAL 


Fic. Energy distribution 566 fission fragments from uranium 233 in. 
The depth the chamber was equivalent 5.3 mm. path standard nitrogen. 


662 fragments. 2.3 mm. path nitrogen. 

680 fragments. 1.1 mm. path nitrogen.* 

1274 fragments. 1.1 path nitrogen. 
While the spectra and were obtained the photographic method, 
that was obtained with electronic pulse analyser. the was 
not calibrated Mev. the are arbitrary 


spectra the energy lost the two groups fragments cannot resolved. 
definite separation into two groups can made only greater equivalent 
depths the chamber. For instance, Fig. shows the spectrum pulse 
heights the first 8.8 mm. track. Here the two groups give energy 
losses that are sufficiently different for resolution. 


*Two experimental points have unfortunately been omitted. These lie just outside the upper 
border the graph. 


100 
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The spectrum pulses obtained with depth chamber equivalent 
1.1 mm. nitrogen shows interesting asymmetry (Fig. C). 
similar asymmetry also evident for depth equivalent 2.3 mm. 
nitrogen (Fig. The spectrum rises more steeply the low energy side 
contrast the behavior depths 5.3 mm. (Fig. and 8.8 mm. 
(Fig. B). This point was confirmed repeating the measurement 
the spectrum for depth equivalent 1.1 mm. using electronic pulse 
analyzer (7), which was capable giving better resolution than the photo- 
graphic method. The spectrum pulses shown Fig. and its shape 
confirms the existence the features observed using the photographic method. 
not possible account for this behavior terms source thickness 
energy losses would result more gradual rise the low energy side 
the peak. Consequently the observed spectra indicate the following. 

The low energy group fragments, represented the left hand peak 
the spectrum total energies shown Fig. gives rise broader 
spectrum than does the high energy group. The ratio the most probable 
energy the low energy group that the high energy group about 
0.65, agreement with earlier work (9). 

Fig. (equivalent depth chamber 8.8 mm. nitrogen), 
the two peaks are closer together. the left hand peak again identified 
the low energy group, the spectrum representing the low energy group still 
broader representing the high energy group. The ratio the 
energies spent the two groups now greater and equal about 

Fig. (equivalent depth chamber 5.3 mm. nitrogen), the 
two peaks are longer resolved, but the more gradual rise the left hand 
side the spectrum indicates that the ratio the energies spent the two 
groups still less than unity. 

2.3 mm. and 1.1 mm. nitrogen), the two peaks are still unresolved, but the 
more gradual fall the right hand side the spectrum indicates that the 
ratio the energies spent the two groups fragments has become greater 
than unity. other words, the ‘unresolved peaks’ representing the two 
groups, which have gradually become closer together with decreasing pressure, 
have crossed over. physical terms this means that the initial ionization 
the low energy group greater than that the high energy group. 

Finally, Fig. plotted the mean energy spent ionization fission 
fragments various distances from the beginning the track obtained from 
these distributions. The mean energy spent the first millimeter the 
track nitrogen N.T.P. 6.5 Mev. 


Discussion 


The conclusion that the initial ionization greater for the low energy group 
fission fragments than for the high energy group accord with the results 


j 
} q 
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measurements the specific ionization fission fragments made 
Lassen (14), which were published after the completion the present work. 
There disagreement this point among the various theoretical estimates 
the initial ionization that have been made. The ratio the energy spent 
the first millimeter track the low energy group fragments that 


MEAN ENERGY SPENT (MEV) 


NITROGEN) 


Fic. Graph mean energy spent the fragments (both groups) function 
the distance from the source nitrogen 


spent the high energy group according Bohr (2) less than unity. Lamb 
(11) and Knipp and Teller (10), the other hand, obtain values for the ratio 
that are greater than unity. The results the present work well those 
Lassen mentioned above are thus qualitative agreement with the predic- 
tions Lamb and Knipp and Teller far the ratio initial ionization 
the two groups should mentioned that evidence conflicting 
with the results Lassen and the present work has recently been published 
Demers (6) and Sherr and Peterson (17), who obtain higher initial 
ionization the high energy group fragments. 
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SURFACE AREA MEASUREMENTS FINE POWDERS USING 
MODIFIED PERMEABILITY EQUATIONS! 


Abstract 


The several modified Kozeny equations for the measurement the specific 
surface fine powders gas-permeability methods have been compared. The 
equations all give very similar surface area values for set experimental data. 
The physical significance the empirical constants used discussed and 
preference indicated. 


Introduction 


Within the past two years, two very similar equations have been proposed 
modifications the Kozeny equation for the measurement the surface 
area very fine powders gas-permeability methods. The use the 
Kozeny equation for the determination the specific surface powders 
the permeability method was first proposed Carman (5) and was found 
give values for the specific surface powder for particles having specific 
surfaces less than about 5,000 sq. cm. per cc. that were satisfactory 
agreement with other methods, such sedimentation and microscope counting. 
was found, however, that, the particle size decreased and the specific 
surface increased, divergences were observed between the gas-permeability 
value specific surface and that obtained other methods. These diver- 
gences have been accounted for modifying the Kozeny equation account 
for the ‘slip’ that occurs when gas flows through capillaries the same 
order magnitude the mean free path the flowing gas. 

comparison two similar modifications the Kozeny equation pre- 
sented this paper. 

Theoretical 


The Poiseuille type equation for capillaries any shape takes the form (4), 


where actual velocity the fluid through the bed centimeters per 
second, 
mean hydraulic radius centimeters, 
pressure drop across the bed dynes per square centimeter, 


shape factor, whose value generally lies between and (4), 

viscosity the fluid poises, and 

actual length path traversed the fluid through the bed 
centimeters. 


received November 26, 1947. 
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For circular capillaries, r/2, where radius the capillary, and 
and Equation (1) takes the form the Poiseuille equation. 
Application this equation granular beds (4) putting 


leads form the original Kozeny equation (7), 


(2) 
where approach (apparent) velocity the fluid centimeters per 
second, 
volume pore space per unit volume (porosity), 
length bed centimeters, 
rate flow the fluid cubic centimeters per second, 
cross-sectional area the bed square centimeters, 
specific surface the granules square centimeters per cubic 
centimeter, and 
Kozeny’s constant, which approximately equal 
Rigden (9) applied the same considerations used above the classical 
Poiseuille equation corrected for ‘slip’, and arrived equation having the 


e 


where having the value 0.874, and 
mean free path the gas molecules. 


Carman (6) has since pointed out that the slip correction factor used 
Rigden was the form used for circular capillaries, whereas, for granular 
beds, the corréction should larger than that used Rigden. recom- 
mends the use slightly modified equation used Lea and Nurse (8, pp. 
27-33), which may written: 


the fraction molecules undergoing elastic collisions with the 


surface, and 
(pi mean pressure the gas dynes per square 
centimeter. 


{ 
| 
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Arnell (1), analogy with the Knudsen semiempirical equation for the 
flow gas through capillary tubes all pressures, arrived the following 


equation: 


where gas constant absolute units, 
absolute temperature, 
molecular weight the gas, 
variable factor, having value approximately 0.9. 
fraction effective void area and given the equation, 
1.40. 


From comparison Equations (4) and (5), seen that the first terms 
the right-hand sides are equal and identical the simple Kozeny equation 
(Equation (2) The second terms the right-hand sides should also 


equal, hence: 


which reduces to: 

€ t),/% Roe = koa. 


taken equal 0.84 (6) and 0.9 (1), 


and actual equality obtained, 0.79, which quite reasonable value. 
Thus, Equations (4) and (5) are identical, 


(6) 
Arnell (1) has defined function which could expressed the 
equation 1.40, and was therefore variable, which 


increased with increasing porosity and approached the limit. Carman (4) 
defined ko/k equal that is, attributed this factor the 
tortuous path that must followed fluid traversing granular bed. 
Within the normal working range porosities, has been found have 
value that lies between 0.4 and 0.6, and ko/k should constant for any 
bed where the ratio the actual length path traversed, the length 
the bed, isa constant. now necessary re-examine available data 
determine whether (or its equivalent) increases with increasing porosity 
remains approximately constant. 


Calculations 


possible study this aspect the permeability equation recal- 
culating the data quoted Arnell (2) for the flow air through beds finely 
powdered quartz and inorganic pigment powders. originally evaluated 
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using his empirical relation between and and calculated the specific 
surface using only the second term the right-hand side Equation (5). 
The specific surface can also calculated using the Kozeny part (first term) 
either Equation (4) and (5) and then applying the result give value 
This has been done here and the results are given Table together 
with the original specific surface values. 


TABLE 


VALUES THE SPECIFIC SURFACE CALCULATED DIFFERENT METHODS FROM THE 
EXPERIMENTAL DATA OBTAINED ARNELL (2) BEDS FINE POWDERS 


Specific surface Specific surface 
(original), recalculated), 
Powdered quartz 
0.425 9.4 0.40 
0.428 9.3 0.43 
0.439 8.5 0.48 
0.475 8.0 0.50 
0.475 0.40 
0.498 10.3 0.38 
0.510 10.1 0.39 
0.514 9.4 0.40 
0.517 9.7 0.40 
0.528 10.0 0.40 
0.537 8.7 0.46 
0.545 9.0 0.44 
0.563 7.8 0.52 
0.584 8.9 0.42 
0.585 9.4 0.39 
0.604 8.8 0.41 
0.610 9.1 0.45 
0.619 0.43 
0.621 7.6 0.51 
0.43 
Mapico Red No. 297 
0.636 46.3 0.61 
0.677 0.63 
0.692 43.2 0.54 
0.718 0.56 
0.721 40.2 0.55 
0.752 0.55 
0.766 39.8 0.61 
38.7 0.58 
Mapico Yellow Lemon 
0.676 29.6 0.48 
0.700 28.4 0.51 
0.704 25.8 0.60 
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TABLE I—Concluded 


VALUES THE SPECIFIC SURFACE CALCULATED DIFFERENT METHODS FROM THE 
EXPERIMENTAL DATA OBTAINED ARNELL (2) BEDS FINE POWDERS—Concluded 


Specific surface Specific surface 
(original), (recalculated), 
Mapico Yellow Lemon—Concluded 
0.754 24.1 0.63 
0.767 0.56 
0.780 0.62 
0.807 21.4 0.74 
26.4 26.3 0.59 
Hydrated yellow iron oxide 
0.702 33.9 0.56 
0.721 33.9 0.52 
0.747 29.3 0.60 
0.771 31.6 0.54 
0.778 28.5 0.53 
0.781 27.9 0.64 
30.9 0.57 
Hydrated chromic oxide 
0.656 
0.664 
0.682 
0.691 
0.706 
0.711 
Pure red iron oxide 
0.479 12.4 0.48 
0.543 12.6 0.42 
0.552 0.44 
0.568 10.5 0.51 
0.583 11.0 0.48 
0.45 
Pure chromic oxide 
0.464 14.9 0.34 
0.487 0.35 
0.493 13.8 0.39 
0.521 7.7 0.43 
0.526 0.40 
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Discussion 


Lea and Nurse (8, pp. 27-33), discussing the modified equations, expressed 
physical picture, but used slip factor, which comparison with the 
equations previously discussed, gives the following equivalents: 


Equation (5) 7.65 (7) 

Lea and Nurse have calculated the values the slip factor, for several 
powders. interest consider these values and convert them 
values means Equation (7). These valuations are reproduced 
Table II. 


Equation (4) 


TABLE 
SLIP FACTORS AND VALUES FOR VARIOUS SUBSTANCES 
(From Lea and (8, pp. 27-33)) 


Substance Porosity Slip factor, 

Cement 240 0.450 3.29 0.43 
0.475 0.43 

0.500 0.46 

0.525 0.43 

Cement 299 0.475 0.43 
Zinc dust 0.425 3.63 0.47 
Lithopone 0.500 4.07 0.53 
0.540 3.53 0.46 

0.600 2.94 0.38 

0.650 3.77 0.49 

Kaolin 0.500 0.36 
0.600 2.78 0.36 

0.625 0.42 

0.650 2.91 0.38 


Rigden (9) found that his equation (Equation (2) did not give constant 
value for the specific surface powder over range mean pressures, 
but showed fall-off the value for the surface low pressures. This effect 
due the use constant, having too value. reading from 
his curves Figs. and and using Table his paper, possible 
recalculate Rigden’s data using either the Lea and Nurse type equation 
Equation (5) and thus obtain constant values for the specific surfaces. The 
resulting values are given Table III. 

From the calculations Tables II, and III, appears that the values 
lie, general, within the range 0.4 and 0.6 and show real trend 
with increasing porosity. This the same range was found for ko/k 
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TABLE III 
SLIP FACTORS AND VALUES @/€ FOR VARIOUS SUBSTANCES 
(After Rigden (9) 


fic Specific surface 
Portland cement 1.28 4.3 0.56 1.41 
Rapid-hardening cement 3.4 0.44 1.58 
Fine slate powder 2.80 3.9 0.51 3.05 
‘Air-blown’ slate 8.43 0.46 9.44 


and therefore seems reasonable assume that Equation (6) true 
equality. The term expressed terms the increased length the 
tortuous path, followed fluid traversing bed granules, over the 
actual length the bed, and, appeared the derivation the original 
Kozeny equation, should appear the modifying term. Also the ratio 
does not show definite increase with increasing porosity, should 
definition, and therefore appears advisable substitute for a/e 
Equation (5) and use the following form: 


where has probable value about 2.5. 

now interest re-examine the permeability data carbon blacks 
reported Arnell (3) and compare the specific surface values these blacks 
given Equation (8) with those quoted his paper. The carbon blacks are 
particular interest that the majority the beds examined showed 
constant mass flow over fairly wide mean pressure range, that only the 
second term the right-hand side Equation (8) enters into the calculation 
specific surface. 

Table are presented the recalculated values the specific surfaces 
the carbon blacks using 2.5 and 0.96 (3) Equation (8), together 
with the original values obtained from the permeability data and electron 
microscope counting. Where the first term the right-hand side Equation 
(8) equation) could applied, the data are included. 

general, the specific surfaces the carbon blacks calculated means 
Equation (8) agree very closely with the original calculations. This fact 
lends support the view that probably equivalent ko/k and may 
replaced it. Further experimentation should carried out settle 
this point, but for the present either Equation (4) Equation (8) may 
used for the calculation the specific surface very fine powders where 
‘slip’ occurs the gas flow measurements. 
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TABLE 
SURFACE AREAS CARBON BLACKS 
(From Arnell (3) 


Specific surface, sq.m./cc. 


Black Porosity Permeability 
Electron 
Original Equation (8) Kozeny microscope 
0.579 
0.638 
0.610 105 
0.636 
0.645 102 
Kosmobile 0.638 155 156 
0.670 150 143 
0.672 155 146 148 141 161 
0.673 131 125 
0.674 131 125 
0.675 155 147 
0.715 
0.721 
Dixiedensed 0.696 164 150 118 
0.704 144 150 129 134 165 
0.716 141 124 
Kosmobile 0.653 178 175 
0.666 197 182 189 176 162 
0.673 171 163 
0.692 126 116 
0.718 
Spheron 0.676 189 179 
0.696 180 184 164 170 144 
0.696 182 166 
Voltex 0.707 344 307 
0.717 328 326 288 280 220 
0.765 306 245 
Neo Spectra Mk. 0..823 435 309 
0.825 418 428 296 303 259 
0.828 430 303 
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SOLAR ECLIPSE OBSERVATIONS THE IONOSPHERE! 


Abstract 


The results are presented ionosphere measurements made the region 
totality during the eclipse July 1945. analysis the results shows 
effective recombination coefficients 1.6 for region, 1.4 for 
region, and 1.0 for region. Evidence variable coefficient for 
region agrees with the suggestion Massey and others high negative ion 
density this level. 


Introduction 


Because the presence the ionized regions the atmosphere, commonly 
known the ionosphere, due mainly absorbed ultraviolet energy from 
the sun, solar eclipse affords unusual opportunity study the effect 
known change radiation over short period. 


For the 1945 total eclipse Canada, site was chosen Victoria Beach, 
Manitoba, the shores Lake Winnipeg. Fig. general view showing 
the lorry, which housed the equipment, the antenna mast, darkroom hut, and 
supplies tent. This point (lat. 50°48’N, long. 96°32’W) was directly beneath 
the totality belt the 200 km. level. Obscuration the sun here was about 
95% the ground, 98% 100 km., and 100% 200 km. The eclipse 
occurred the early morning, the sun’s elevation being less than first 
contact and about 25° fourth contact. 


Observations 


automatic ionosphere recorder was constructed the National Research 
Laboratories especially for the eclipse observations, and the equipment was 
set the site and operation July Details this equipment are 
found the N.R.C. Radio Branch Report PRA No. 2isa 
front view showing the rack the left, the scanner and 
visual monitor center, and the camera recording unit the right. The equip- 
ment produced virtual height versus frequency curves mm. film, 
covering the range Mc. per sec. Mc. per sec. logarithmically just 
over two minutes. Following observational program laid down the 
Canadian Radio Wave Propagation Committee for the ionosphere stations 
Canada, records were obtained min. intervals throughout the month 
from July Aug. except the eclipse day when records were made 
five minute intervals over eight hour period (Fig. 3). 

Conditions were quite disturbed for few days preceding the eclipse when 
strong auroral displays were seen. Data for the eclipse day are compared 
with control data for three days before and three days after the event. Owing 


Manuscript received October 1947. 
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slight disturbances the region July and 12, data for these two 
days were eliminated region calculations. 

The photographic records for the days July inclusive have been care- 
fully analyzed. Virtual heights the ionized regions were read accuracy 
km. and critical penetration frequencies +0.05 Mc. per sec. The 
actual time observation any particular critical frequency was tabulated 
the nearest half-minute. 


Analysis Results 


The object the following analysis investigate the general behavior 
maximum electron density each the regions F,, and The 
correlation between the observed values maximum density during the 


eclipse and the values predicted the simple recombination law, 


studied for the three regions. Maximum electron density, 
derived from the equation relating the critical penetrating frequency 
the ordinary ray, 1.24 where maximum electron density 
the region terms electrons per cubic centimeter, ionizing rate the 
level maximum density terms electrons per cubic centimeter per 
second, recombination coefficient, and frequency cycles per sec. 

The following the method analysis used interpreting the eclipse 
effects the ionosphere. From the control data normal day electron 
density curve for July was drawn interpolation for each the three 
recombination coefficient, for each the three regions the curve ionizing 
rate, was drawn Curve Figs. Making further assump- 
tion that the ionizing rate given time during the eclipse was proportional 
the visible area the sun’s disc, Curve was constructed for each region. 
Finally using values from Curve the initial value and from 
Curve first contact, and the assumed value curves were constructed 
values were assumed and curves constructed although only those lying 
close the observed points are shown. 

Let first consider the correlation shown for region results Fig. 
Because does not regularly appear until around a.m. July, Curve 
extrapolated beyond the mean the observed values, back the time 
the beginning the eclipse. Curves and are electron density curves cal- 
culated from assumed values 0.7 and 1.4 Observed 
points for maximum electron density are fairly close agreement with the 
curve for 1.4 except near the beginning the eclipse. The 
deviation here most probably due partly error extrapolation the mean 
data and partly the normal day day spread can seen 
that the minimum value for and the time after totality which 
occurs correspond quite closely the curve 1.4 
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region results not show the same degree correlation observed 
for Fig. shows derived curves for two assumed values 1.6 
and 1.9 Curve smooth curve showing the general trend 
observed points. Although Curve has the same minimum falling 
about the same delay time the curve 1.6 interpretation 
this limited the accuracy the observations qualitative one, but 
the discrepancy may either the assumption that the sun’s radiation 
proportional exposed area only, the assumption that the region obeys 
simple recombination law. one assumes that the sun’s corona contri- 
butes considerable proportion the ionizing radiation cause more 
rapid fall and rise near first and fourth contacts, then one would expect 
considerable ionizing rate totality from the annular ring the corona which 
still exposed. allow electron density fall low that observed 
experimentally large value must then assumed. This, however, 
would not result the same order delay time for the minimum value 
that observed min.). This alone would seem rule out any large 
coronal effect, apart from the fact that none definitely observed results 
the other regions. 

the interpretation that the region obeys different recombination 
law from the simple one assumed, measurements the observed curve show 
large change effective recombination coefficient from near the time 
first contact totality. The measurements fit well with the conclusions 
Massey, Bates, and others (1, who show that there probably high con- 
centration negative oxygen ions region levels. This results 
effective recombination coefficient dependent the ratio ion density 
electron density neutral medium (2). The equation equilibrium 

iven (1) 

where electron density, electrons per cubic centimeter, 

positive ion density, ions per cubic centimeter, 

ratio negative ion electron density, 

recombination coefficient electrons and positive ions, and 

the coefficient for electron transfer from negative positive 

ions. 

which are derived from controlled laboratory experiments. Because 
large may simplify Equation (1) obtain the approximation 
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the time 1127 1/2 G.M.T. the value reaches maximum 


electron ratio this time. Measurements Curve Fig. give values 
this time 4.50 electrons per cc., 3.58 electrons per cc. 


electrons per cc. per sec? 


2 


Thus, 


(the effective coefficient), 


4.2 


The figure may correct only order magnitude but serves 
useful comparison with the result obtained totality (1215 G.M.T.), 


2100. 


when, being zero, where that occurring 1215 G.M.T. 


From the observed curve, 3.1 electrons per cc. and 


—7.69 electrons per cc. per sec. Thus 


The drop from 2100 400 indicated here shows that the term 
+ 


region results are shown Fig. the eclipse day and the control 
days displayed much more irregularity critical frequency than did 
either the other two regions. The experimental results are, therefore, 
much less significant this case. Curve the mean control day density, 
shows hump about 1230 G.M.T. the records for all the 
days the eclipse month showed only slight possibility systematic 
rise critical frequency this time. The apparent rise fact may 
due only the difficulties analysis, appears shortly after the separation 
regions. was, therefore, assumed that the hump was random 
variation effect, not smoothed out control day data, which could not 
accepted characteristic the region. Curves and were derived 
from Curve (neglecting the hump) for assumed effective recombination 
coefficient 0.1 can seen that the observed points, though 
badly scattered, show some correlation with the curve. Electron densities 
were lower both before and after the eclipse than those predicted from the 
control days, but are well within the expected day day spread. 


oot! 


3$41793 YOS ALIGN3IO O3AYISBO 

AWO 3541793 WOS Z ONIZINO! € 

24 364 WI $4939493 3541793 


OnIZINO! 


= 
2 
~ 
Q 
=) 
i=) 
< 
S 


40 


3961793 40 


3541793 40 OND 


3 
. 
. 
= 
. 
a 
. 
« 
4 
= 


144 CANADIAN JOURNAL RESEARCH. VOL. 26, SEC. 


interesting point noted the absence reflections for the 
time between 1251 G.M.T. and 1321 G.M.T. when maximum densities 
recovered rapidly enough approach the densities this interval. 
Absorption was too great allow observation echoes during this 
time. 

Conclusions 

From the results derived from the eclipse observations Victoria Beach 
some conclusions regarding the general behavior the three main regions 
the ionosphere may drawn. 

The region behaves though has constant recombination coefficient 
approximately 1.4 seems unlikely that negative gas ions are 
formed stably the height 200 250 km. 

The region apparently has varying effective recombination coefficient, 
the average value during the eclipse being about 1.6 making 
some approximations, the value coefficient arrived about 4.2 
soon after first contact, and 0.8 totality. This shows drop 
ratio negative oxygen ions electrons from 2100 400 totality. 

The region shows distinct eclipse effect but the quantitative analysis 
limited the experimental evidence. effective recombination coefficient 
the order seems hold for this region. 

may said that there experimental evidence coronal effect 
the solar eclipse, and therefore the ionization due the corona must small 
compared with that due the main body the sun. 
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HIGH TEMPERATURE CONTROL 


Lutz AND JOHN Woop? 


Abstract 


This electronic circuit useful the control electric air-bath furnaces 
having long time lag. The control accomplished means relay that 
shunts portion the resistance series with the furnace winding. The ther- 
mometric device thermocouple connected series with mirror-type galvano- 
meter and variable voltage supplied potentiometer. The setting the 
potentiometer selects the temperature which the device will control. light 
beam falling the mirror reflected one two photocells. One photo- 
cell closes the shunt means the described circuit and increases the current 
the furnace. This condition maintained until the beam moves across the 
second photocell, when the shunt again opened. The device will control 
temperatures over 1000° within precision over any desired length 
time. 


Introduction 


The temperature control described here was designed for use over 
wide temperature range with air-bath furnaces, electrically heated covered 
resistance elements. long time lag inherent property this type 
furnace. There are number temperature controls described the liter- 
ature (1, which will give higher degree precision but have limitations 
temperature range, power input the furnace, time lag. 

The limitation temperature range usually imposed the range 
the platinum resistance thermometer. When thermocouple used the 
thermometric device, the limit determined the refractory nature the 
elements employed. Thus for tungsten—molybdenum couple, temperatures 
excess 2000° may controlled. 

The furnace resistance wire also determining factor the upper tem- 
perature limit. The furnace present being thermostatically controlled 
has nichrome element permitting operation 1000°C. However, 
furnace constructed employing tungsten windings, thus permitting 
temperatures excess 2000° 

the previous circuits, thyratron was used series with the furnace 
current regulator controlling the relative phase the grid voltage. 
The use thyratron has the advantage continuous control. However, 
such circuits are not applicable over wide range power consumption nor 
they lend themselves readily either frequent changes the required 
temperature changes the furnace windings. Also, the thyratron 
circuits depend upon the light intensity falling single photocell. felt 

difficulties maintaining constant light source over long periods 
time might arise and that, with single photocell, there danger that the 
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light will move completely off the photocell there any appreciable time 
lag the furnace. the latter occurred all control would lost. 

The circuit described uses relay shunt portion (Fig. 
the series resistance and thus affords practical limitation the power 
delivered the furnace. However, gain this respect continuous control 


FURNACE 


Fic. 


has been sacrificed; the effect discontinuous control may minimized 
keeping the shunted portion the series resistance small possible. The 
selection new controlled temperature has been simplified two manual 
adjustments; the setting the series resistance, for approximate tem- 
perature and the adjustment the opposing e.m.f. for the thermocouple 
(Resistance R). 

The chief advantage this circuit use two photocells, one being em- 
ployed open and the other close the shunt. the light beam passes 
either cell, the proper action accomplished does and the light must 
eventually return the opposite cell, ill effects occurring repasses 
the first cell. 

Description the Circuit 


The circuit Fig. consists the following stages: the relay tube, 6F6, 
circuit (2) consisting two 6C5’s, and two pulse circuits 
functioning alternately. Each pulse circuit consists 884 used saw- 
tooth generator, used bias amplifier, and 927 phototube. 

The relay, Aminco ‘supersensitive’ mercury relay, activated ma. 
and releases when the current reduced below current con- 
trolled the bias the 6F6 since this tube series with the relay. The 
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bias may either two values determined the condition the 
Jordon circuit. This circuit (2) has two stable conditions, either which 
one tube conducting while the other cut off. When the right hand 


100 AQ. R; = 50 KQ. Rs = R; = 400 KQ. Rs = 20 KQ. Ry = 25 KQ. Rio = Ru = 
100 KQ. Ri = Rss; = 5 MQ. Ris = Ry = 15 KQ. Ris = Rx» = 25 KQ. Ris = Rox = 


5 KQ. Rie = Roy = 25 KQ. Ri = 20 KQ. Ris = Rio = 100 KQ. Rao = 10 KQ. 


cut off, the bias the 6F6 allows the tube conduct the necessary ma. 
close the relay. When this conducting, the relay current reduced 
below ma. The VR150/30 the voltage divider permits maximum 
change grid bias the same time reducing the drop across the grid 
return resistor. 

The Eccles—Jordon circuit altered from one stable condition the other 
applying negative pulses the grid the tube cut off. These 
pulses are formed passing saw-tooth wave from standard type 884 
sweep circuit through differentiating circuit. Fig. the upper 884 
thyratron supplies negative pulses through the differentiating circuit, 
and Rig, the grid the right 6C5, cutting off and closing the relay. 
similar manner the lower 884 functions open the relay. 

The saw-tooth wave produced only when the bias the thyratron 
lifted point where the tube can bias controlled bridge 
circuit with 6F5 forming variable resistance element one arm. When 


RELAY 
Rio 
927 
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the light falls the corresponding photocell, the negative voltage drop across 
the grid return resistor the 6F5 increased, causing increase the 
resistance this arm the bridge, thus permitting the 884 function. The 
6F5 high-gain tube and care must taken prevent oscillation the 
leads from the photocells are long. Since this d-c. amplifier, permissible 
use large by-pass condenser from plate ground. 


Operation the Circuit 


high degree precision not claimed for this circuit. was designed 
for purpose requiring precision only +2°C. high temperatures 
(900° 1500° C.), which accomplished. necessary this could improved 
considerably increasing the current sensitivity the galvanometer and 
improving the optical system. The device has seen 1000 hr. almost con- 
tinuous operation and has given excellent performance throughout the period. 
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THE PHOTOELECTRIC TEMPERATURES THE CYGNI 


Abstract 


The photoelectric temperatures Cygni stars spectral type from 
have been determined from emission lines Zanstra’s method. 
approximately half the stars, accurate spectrophotometric measures emission 
line intensities were available and for the remainder the intensities were estimated 
the basis published descriptions spectra. The computed temperatures 
were found range from 10,000° 42,000° and after smoothing graphical 
methods temperature scale for the and type stars showed the following 
values: AO, 16,000°; BO, 40,000° Comparisons are made with 
temperatures derived other methods and reasons are given for attributing 
greater weight the values from than those the lower range 
the sequence. 


Introduction 


There are three principal methods now use determine the temperatures 
stars, (A) the application Planck’s formula the observed distribution 
energy the continuous spectrum (6), (B) the use ionization formulae 
such those Saha conjunction with observations relative line inten- 
sities atoms different stages ionization (10), (C) the determination 
photoelectric temperatures the basis the intensity emission lines relative 
the continuous spectrum fourth method makes use the Boltzmann 
formula for the distribution atoms excited states but limited 
application among the early type stars (13). the present paper data 
the intensities emission lines have been used derive photoelectric 
temperatures for number stars Cygni type. 


Theory Method Used 


Method (C), developed Zanstra (14), was originally used deter- 
mine the temperatures planetary has since been adapted 
include stars that have around them envelopes such may have their 
origin the ejection atoms from the surface the star. This has been 
done for the Wolf Rayet stars, for number stars Cygni type, and for 
several novae Beals (1, 3), for Nova Lacertae McKellar (7), 
have this method available the case the Cygni stars since Method (A) 
invalidated the distortion the continuous spectrum due space 
absorption, and Method (B) made difficult application the combination 
emission and absorption Cygni type spectra. 
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The application Zanstra’s method involves the following assumptions, 

(1) The star assumed black body radiator temperature 
obeying Planck’s radiation law. 

(2) The envelope the star assumed consist atoms one element 
only, and majority the atoms are assumed the normal state. 

(3) The ultraviolet starlight emitted the star completely absorbed 
beyond the ionizing frequency the atoms present the envelope, causing 
the ejection photoelectrons. 

(4) The photoelectrons recombine with the ionized atoms the upper states 
the atom and steady state the number quanta ultraviolet starlight 
absorbed equal the number recombinations any time interval. 

The observational data required for the calculation are values the 


quantity where 
A = Le (1) 


where represents the total energy the emission line, average frequency 
and the total energy per frequency unit the continuous spectrum for 


the same frequency. has been shown Beals (1, that this quantity 
(2) 
where the are the velocity displacements from the center the band 
the edges, the intensity the band relative the continuous spectrum, 
and the velocity light kilometers per second. The quantity 
the equivalent width the line expressed units velocity displacement. 
the basis the above assumptions the theory given by-Zanstra leads 
the equation 


— 1 dx = A. (3) 
where 


constant, the wave length angstroms, the frequency, 
the ionization potential volts, and the ionizing frequency and wave 
length respectively, and the summation carried out over all the visible 
emission lines for particular atom. The equation solved trial order 
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determine the value The temperature determined the ‘photo- 
electric’ temperature referred the title this paper. more detailed 
discussion the theory and its application individual observations has 
been given Beals (3). 


Observational Data and the Use Estimates Intensity 


For considerable number the stars considered, spectrophotometric 
measures equivalent width made Victoria were available. These values 
equivalent width were derived from plates calibrated with the aid 
rapidly rotating step sector (100 revolutions per second) placed front 
the slit calibrating spectrograph. The plates were analyzed with the 
aid microphotometer tracings and the process reducing the tracings 
equivalent widths was greatly facilitated new intensitometer, which 
has been described elsewhere The actual observations will discussed 
more detail forthcoming publication the Dominion Astrophysical 
Observatory. illustrate the general type spectrum dealt with this 
paper, portions the spectra two Cygni stars have been reproduced 
Plate along with that ordinary absorption line star early type for 
comparison. 

addition those stars for which quantitative measurements line 
intensity are available there are numerous others for which there are only 
general descriptions and some cases published reproductions spectra. 
Ordinarily might supposed that such instances significant calculation 
temperatures could made. However, Zanstra has previously pointed 
out, this method determining temperatures very insensitive observa- 
tional error. This point may illustrated the following examples;— 

(1) For the star H.D. 31293 taking the equivalent width the line 
hydrogen 200 and using standard values for the ratios 
and the temperature the star was found Then 
taking the value the equivalent width 400, twice much before, 
and using the same ratios the value 13,500° was found. 

(2) For higher values the temperature the effect more noticeable but 
again not very the star H.D. 151804, using the measured value 
the equivalent width and doubling the value, the tempera- 
tures 40,800° and 42,800° were calculated, whereas, for the same star, 
measures different line from different ion gave temperatures 32,400° 
and 34,600° Since the method itself and the assumptions which 
based involve uncertainties the order few thousand degrees for individ- 
ual stars, the errors due the use estimates are probably not very 
significant. Accordingly estimate line intensity has been made for every 
for which available and believed that reasonable 
degree confidence may placed the temperatures calculated from these 
estimated intensities. 

our calculations the temperature was calculated for both measured 
judged and for twice A,. This the same assuming first, that the 


CS 
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number observed emission quanta equal the ultraviolet quanta 
absorbed, and second, that the observed are only half those absorbed. The 
first temperature may regarded minimum temperature, and will 
denoted min. our results. 


should stated here that the temperatures obtained for the same star 
using different atoms ions show large differences, result that naturally 
leads doubts the accuracy the method. These differences will 
discussed later part the paper. 


Evaluation the Integrals 


For use our calculations require tables the integral 


x? 
qi = 7-1 dx 


Xo 
Zanstra (14) gives table this integral for values between and 
11.0 calculated Debye’s method. our work the value some- 
times reaches values over that extension necessary. This can 


where the ionization potential volts and the absolute temperature 
15,000° less reach values near 20. 

calculating values for the integral beyond equal approximation 
was made that will shown justified. require 

oO 

Xo 


xo 
Integrating parts: 


testing this approximation obtain the values Table 


TABLE 
Xo [ ve xe 
To 
0.00554 0.00554 


0.2493 0.2500 


[ 
Xo 


PLATE 


Call 
Call Hel SilV 
PLATE 


Emission Line Stars Cygni Type. 
(3) Ordinary absorption Line Star For comparison. 


— - - - — — 
Call Fell 


”™ 


BEALS AND HATCHER: PHOTOELECTRIC TEMPERATURES THE 


thus see that even for values quite low the table the approximation gives 
very good results. 
Similarly may calculate values the second integral used Zanstra. 
This 
x3 


Zo 


Zo 


which becomes 


using similar approximation. 
Upon integration parts before this gives 


Trying obtain 1.590 this method, whereas the value 
Debye’s method 1.594. 

shall however need numerical values for the first integral only and these 
calculated according the formula above are shown Table II. 


TABLE 
Zo Zo 

10.0 0.00554 16.0 0.0000326 
11.0 0.00242 17.0 0.0000135 
13.0 0.000445 19.0 
14.0 0.000188 20.0 
15.0 


may noted that Mohler (9) gives table calculated expansion 
method for values 15.5 agreement with those given here. 


Use Graphs 


facilitate computation, graphs were drawn from which several the 
values appearing Equation (3) could read off directly. 


From (3) 1.432 108 


where the wave length angstroms which for the range dealt with here 
varies from for any definite value 


may calculated giving values 


between and This was done initially for values 10, 15, 
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20, 25, 30, 35, and thousand degrees and this set seven curves was used 
the first few calculations, being considered sufficient interpolate 
between these obtain the temperature sufficient accuracy. Later, 
however, when was noticed that most temperatures lay the 10,000° 
20,000° group, extra curves were drawn for Tequal 12,500° and 17,500° 
and interpolations made between these values. Occasionally greater accuracy 
was used for special stars and interpolations between 1000 degree intervals 
were made. such cases the values obtained were checked with those 
obtained interpolating over wider intervals and was noted that they 
agreed sufficient degree accuracy. cannot claim therefore that 
the values given here are precise could obtained drawing many 
more intervening curves but felt that the accuracy amply sufficient for 


the purpose. 
Treatment Hydrogen Lines 


hydrogen played important role the spectra these stars method 
shortening various calculations was devised, using assumed value for 
the Balmer decrement. 

will noted that the equivalent widths for the main Balmer lines 
the star Cygni are 


leading the ratios Several the other stars 
showed roughly these ratios and where only had been measured 
was assumed that the ratios held, and calculations were made accordingly. 
come cases when the hydrogen lines showed values different from those 
given above, separate calculation was made for the temperature. This 
ratio not the Balmer decrement given some authors, but its connection 
with can seen from References (1, (The values given here are each 
case relative the continuous spectrum the appropriate wave length.) 


Classification 


most cases ordinary absorption line intensity ratios are not 
for the Cygni stars because the blending absorption and emission lines. 
The spectral classes included the second column Table have therefore 
been determined from the general level ionization the stellar spectra. 
Dr. Petrie very kindly undertook classify all stars for which data 
were available during visit this Observatory January, 1947, and most 
the designations Table II] are his. more detailed discussion the 
problems the classification the Cygni stars will appear forth- 
coming publication the Dominion Astrophysical Observatory. the mean- 
time believed that classifications are precise enough for the present 
purpose determining general relation between spectral class and tem- 
perature for the Cygni stars. 


— 


ves 
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Table Results 


The results the calculations are given the following table. Column (1) 
shows the star, (2) its classification, (3) the atom ion encountered its 
spectrum emission, and its ionization potential, (4) the wave length the 
line, (5) the equivalent width kilometers per second (6) min. and (7) 
order distinguish between the values equivalent width that have been 
directly measured and those that have been estimated, the latter values are 
enclosed parentheses. 

will noted that many more lines occur for some stars than for others. 
This either because the emission spectrum more completely developed 
because more complete data are available. some instances where 
spectrum shows marked variations the temperature has been calculated for 
more than one epoch, the case the star 108. 


TABLE 


TEMPERATURES THE STARS 


Star Class Atom LP. °K. °K. 
H.D. 108 13.54 6562.817 295 13,400 15,200 
(1934-38) 213 
(1944-45) 6562.817 10,500 11,100 
4650.16 
4651.35 
(1944-45) 5696.0 35,200 36,400 
4650.16 
(1934-38) NIII 47.2 37,200 39,700 
4634.16 
(1944-45) 36,100 38,000 
(1934-38) Hell 4685.682 40,500 42,300 
M.W.C.419 13.54 4861.332 (200) 15,200 16,100 
Fell 16.16 4351.764 (75) 14,500 15,300 
(75) 
H.D.7902 6562.817 100 10,700 11,400 


_ 
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TABLE 


TEMPERATURES THE STARS—Continued 


H.D.12953 13.54 6562.817 10,300 10,800 

(270) 
650 
5015.675 (98) 
4471.477 (70) 
Fell 16.16 (75) 14,500 15,300 
4233.167 (75) 
H.D.21389 13.54 6562.817 (100) 10,700 11,400 
H.D.31293 13.54 (200) 12,800 13,500 
H.D.31648 13.54 6562.817 (200) 12,800 13,500 
H.D.32256 13.54 6562.817 (224) 12,600 13,200 
(13) 
lanic Cloud (325) 


4101.737 (205) 
3970.074 (129) 


4101.737 
M.W.C.497 Fell 16.16 4351.764 (75) 14,500 
(75) 
H.D.41511 13.54 6562.817 360 13,400 14,600 
H.D.45910 
(1937) 13.54 158 14,700 15,800 
(Av. 1939-44) 6562.817 1730 16,100 17,700 
(25/3/39) 6562.817 815 14,900 16,000 
(1943) 13,500 14,700 
(Av. 1923-24) 16.16 4351.764 88.3 15,000 
(1926) 4351.764 30.1 14,200 
4233.167 30.1 
(1943) 4351.764 13,800 15,100 
46.6 
H.D.51480 13.54 6562.817 2740 16,600 18,200 
345 
103 
4101.737 


3970.074 


— 

— 
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TABLE 


H.D.51480 Fell 16.16 15,000 15,700 
4351.764 
4233.167 
4178.855 
4173.450 
H.D.51585 13.54 6562.817 (100) 10,700 11,400 
H.D.53179 13.54 6562.817 (100) 10,700 11,400 
M.W.C.560 13.54 6562.817 (436) 13,200 14,400 
(49) 
4233.167 
(225) 
H.D.94878 13.54 (783) 17,500 19,300 
H.D.94910 13.54 (783) 17,500 19,300 


Hel 24.48 6678.149 (227) 25,900 27,800 
5875.618 (540) 


650 

5015.675 (197) 

4471.477 (139) 
H.D.138403 4861.332 (783) 17,500 19,300 
H.D.151804 Hell 54.17 4685.682 (56) 40,800 42,800 

4634.16 (84) 
47.67 5696.0 (16) 32,400 34,600 

(13) 
H.D.152408 Hell 54.17 4685.682 (28) 39,000 40,800 

4634.16 (20) 


— 


CANADIAN JOURNAL RESEARCH. 


VOL. 26, SEC. 


TABLE 


TEMPERATURES THE CyGNI 


Star Class Atom I.P. 


Fell 16.16 4351.764 (150) 15,300 16,700 
(150) 
Hel 24.48 6678.149 (170) 25,200 26,900 
(400) 
650 
5015.675 (148) 
4471.477 (105) 
Fell 16.16 4351.764 (68) 14,400 15,200 
(68) 
H.D.163296 6562.817 (360) 13,400 14,600 
H.D.168607 4861.332 (392) 16,200 17,800 
H.D.169226 Hel 24.48 6678.149 (170) 25,200 26,900 
(400) 
650 
5015.675 (148) 
4471.477 (105) 
H.D.169454 6562.817 436 13,200 14,400 
Hel 24.48 5875.618 18,100 18,900 
650 
H.D.169515 Hel 24.48 5875.618 (400) 25,200 26,900 
650 
47.2 (3) 30,000 32,000 
4634.16 (3) 
H.D.174638 Hel 24.48 5875.618 (400) 25,200 26,900 
650 
H.D.183143 6562.817 (11) 10,100 10,300 
(13) 
Hel 24.48 5875.618 17,900 18,600 
650 
4634.16 10.7 
H.D.190073 6562.817 1304 15,400 16,600 
190 


3970.074 


| 


BEALS AND HATCHER: PHOTOELECTRIC TEMPERATURES THE CYGNI 


TABLE 


TEMPERATURES THE STARS—Continued 


H.D.190073 Fell 16.16 5018.434 15,800 
4923.921 
4555.890 
4522.634 
4515.337 
283 
4351.764 
4233.167 
4571.971 12,700 13,300 
12.8 4246.829 10.000 10,200 
H.D.190603 13.54 6562.817 224 12,600 13,200 
650 
Cygni 13.54 6562.817 3660 17,500 
4861 .332 783 
325 
4101.737 205 
3970.074 129 
5875.618 540 
650 
5047 .736 
5015.675 197 
4921.929 
4713.143 
4471.477 139 
4143.759 
4120.812 
3964.727 
3871.819 
3867 
29.49 6482.07 13.4 26,100 27,800 
5747.29 
5710.76 8.5 
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TABLE 


Star Class Atom I.P. 
H.D.193137 
H.D.197345 13.54 
H.D.198478 13.54 
M.W.C.374 
H.D.207757 13.54 
H.D.218393 13.54 
Fell 16.16 
H.D.221650 
Andromedae) Hell 54.17 
47.2 
Companion 
Aquarii Hell 54.17 
NIII 47.2 
H.D.223385 13.54 


E.W., 
5679.56 28.0 
13.8 
19.8 
24.5 
4630.537 32.9 
4621.392 15.0 
14.4 
25.8 
4447 .033 6.6 
18.6 
6.7 
5696.0 1.0 
6562.817 
6562.817 
50.7 
4555.890 19.8 
29.3 
4351.764 29.4 
4233.167 
4178.855 18.5 
4173.450 20.9 
(1566) 
(340) 
5875.618 (810) 
650 
5015.675 (295) 
4471.477 (219) 
6562.817 806 
60.5 
4101.737 
17.9 
123.4 
4640.64 
4634.16 
(80) 
(46) 
4634.16 (20) 
6562.817 (100) 


27,000 
10,100 
9,500 


15,100 


19,400 
26,900 


14,200 


12,500 
43,200 


36,100 


41,100 
36,100 


10,700 


28,700 
10,300 
10,000 


16,700 


15,500 


14,700 


42,900 


11,400 


| 
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Discussion Results 


discussing the validity the results this study important bear 
mind the fundamental limitation the method, namely, that claims 
only set minimum value for temperature star. This limitation 
arises largely consequence two assumptions (a) that all the radiation 
the violet the ionizing frequency atom absorbed the process 
photoelectric ionization and that most electron captures are made 
upper atomic levels. 


Referring assumption (a), even could supposed that all, 
majority the atoms envelope, were one kind might well happen 
the case incipient emission that the envelope would insufficient 
optical thickness for complete absorption and this would lead too low 
estimated temperature for the star. Also mixed atmosphere, owing 
the existence overlapping regions absorption, assumption could 
realized only under very exceptional circumstances, and the failure the 
assumption would again lead too low temperature. 


There thus conflict between the different values temperature derived 
for the same star the use lines different atoms. Since each represents 
minimum value, the highest calculated temperature (ordinarily derived 
from the atom highest ionization potential) should correspond most closely 
the true temperature the star. Also, the numerous instances where 
stars the same spectral type exhibit different emission characteristics and 
therefore different photoelectric temperatures, the differences may attri- 
buted mainly the different optical thickness their respective envelopes. 

The assumption regarding electron capture upper levels also likely 
result too low rather than too high estimates the stellar temperature. 
Here have followed Zanstra assuming that half the actual captures result 
observable emission. view the work Cillié (5) the capture 
coefficients for hydrogen would appear that most captures take place 
the lower levels, the largest number being the ground state. thus 
possible that considerably less than half the actual captures are made 
atomic levels that result emission lines the spectral region actually 
observed and that all the calculated temperatures are somewhat too low. 

Incomplete observation the emission spectrum, which inevitable 
some instances, will have the same result and would therefore appear that 
the results this paper are likely give rather conservative estimates 
the temperatures attained the and type stars. 


the Use Ions other than Hydrogen 


should noted that Zanstra’s theory was developed for the case the 
envelope star containing hydrogen only (3) because the relative ease 
the calculations for the hydrogen atom. have not attempted further 
development the theory for the more complex atoms and consequence 
the relation the observed line intensities the number absorbed 


- 
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ultraviolet quanta more uncertain than for hydrogen. However, the more 
important cases NIII, CIII) have made use single line 
single multiplet, and there would appear little likelihood that use the 
observed intensities would result overestimate the number ultra- 
violet quanta. The original aim the investigation derive minimum 
temperatures therefore preserved, although other considerations arise where 
the mechanism fluorescence involved. 


Fluorescence Effects 


fluorescence mechanism (4) another factor considered 
analyzing these results. For completeness shall give short discussion 
along the lines used McKellar (7) treating the same 
Considering first the cycle for the case which the line 
Hell produced, the diagram Fig. shows quite clearly how this 
line may occur with anomalous intensity, the exclusion other lines 
its series. Because the large intensity \304 many the helium ions will 


121518 


Fic. cycle for and Hell. 


State represented the diagram. Then because close coin- 
cidence with the very intense Ly, line hydrogen the ion may moved 
State the absorption quantum Ly, and falling back will 
produce the visible line while other members the series will absent. 
The Ly, line hydrogen while that for helium \1215.18, 
that seen that the coincidence very close. 

The OIII NIII cycle more complicated and consists two close 
coincidences. Its mechanism may seen from Fig. presence the 
wave lengths NIII unaccompanied the strong lines 
4103 evidence for its occurrence. 


— 
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connection with the above mechanism fluorescence McKellar makes 
the interesting comment that since each quantum \4640 due Bowen’s 
line absorption process follows from one 304 which resulted from the 


LINES 


LINES 
A4097 
AG52 


2p sp 2s 


Fic. Fluorescence cycle for HeII, and 


photoelectric ionization atom, then the quantum should 
considered evidence the presence ultraviolet radiation beyond the 
ionizing frequency rather than NIII. Similarly the \4686 line 
should considered line hydrogen, and considered the calculated 
temperature the star will general much lower. the former case, 
taking band Hell and calculating the temperatures thus, the 
result obtained will not differ greatly from that found for NIII. 


Because the existence the fluorescence mechanism would appear 
that temperatures calculated the basis the lines NIII and Hell 
would somewhat unreliable although means certain that they 
would all cases too high. However since the action the mechanism 
increase the intensity the line may expected that such lines 
will lead temperatures higher than those due solely recombination. 
Fortunately most instances where lines due and NIII appear there 
are also lines due CIII which does not appear involved the fluores- 
cence mechanism. There reasonable agreement between the temperatures 
calculated for CIII and for NIII and Hell although the latter are somewhat 
higher. constructing Fig. was decided make use CIII only 
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plotting the relation connecting temperature with spectral type, and will 
seen that the points for NIII and Hell are mostly above the curve. 


45000 


30000 


TEMPERATUR 


10,000 


SPECTRAL TYPE 
between photoelectric temperature and spectra type for the stars. 


The Stellar Temperature Scale for Early Type Stars 


will noticed that the stars cover the range the spectral sequence for 
spectral types A4. Fig. the photoelectric temperatures found 
have been plotted against spectral type. For the interpretation the graph 
few remarks are necessary. First, the points seen the diagram are 
individual calculations for star for definite ion only and not average 
for several ions cases where several appeared. The curve drawn through 
the highest points each the spectral classes, since the temperature deter- 
minations, has been stated above, refer minimum values, except cases 
where the fluorescence effect may acting. These latter are most noticeable 
the end the sequence and will noticed they lie above the curve. 
Here the NIII and Hell determinations are generally greater than those 
for CIII for which fluorescence effect coincidences line wave lengths 
have been found. The values temperature from NIII and Hell have 
been disregarded drawing the curve. Second, addition the values 
obtained ourselves, there have also been introduced into the diagram 
the values Mohler for the stars, principally and which there 
were few our group. These are indicated open circles. Third, the 
sudden rise temperature the vicinity due principally the 
appearance the spectra measurable lines ions higher ionization 
potential than for example sufficient variety ions was 
available would expect the area under the curve completely filled in. 


40000 
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will also noticed that temperatures derived from hydrogen not show 
much variation with spectral type and the general appearance the curve 
would indicate series steps for the different ions. 

There are several discussions scientific literature the temperature 
scale stars. One the most recent that Kuiper (8), who also gives 
diagram showing the log the effective temperature against stellar sequence. 
Petrie (11) has also given discussion the temperatures the 
stars. The values given for this part the stellar sequence the two writers 
are not agreement, Kuiper giving 25,000° 80,000°, whereas Petrie gives 
28,000° 36,000°. Petrie has pointed out, Kuiper uses values 
obtained Beals (3) for Wolf Rayet emission being applicable the 
associated O-type spectra, whereas now believed that composite spectrum 
really observed and that the emission line temperatures refer the star 
and not the O-type object. 

may seen that for our curve the values obtained are, for BO, approxi- 
mately 30,000°, although this quite rough owing lack stars this 
general region, and, for extrapolated temperature the neighborhood 
40,000°. These figures are slightly higher than those Petrie’s corre- 
sponding classification. For the lower end the scale our values are some- 
what higher than those shown Kuiper. 

There are, well known, some doubts the validity applying 

‘results from such Unusual objects the Cygni stars determination 
the general temperature scale the stellar sequence. However the tem- 
peratures the hotter stars are subject much uncertainty and the results 
different investigators differ widely among themselves that surely 
desirable make full use data from every possible source attempting 
arrive reliable the upper part the scale involving the 
early O-type stars the presence emission lines due Hell, and 
fairly normal feature and shown one the typical stars the group 
H.D. 188001. would therefore seem quite permissible make use these 
lines fix the temperature this end the scale. 

Those early A-type stars that show Cygni characteristics are the other 
hand mostly quite unusual objects. number them have exceptional 
luminosities while one two instances suspected that the spectra are 
composite. consequence the present determination probably entitled 
less weight for the lower end the range than for the O-type stars. From 
number converging lines evidence the temperature scale has been fixed 
10,000° for class and will noted from the diagram that our values 
are approximately 6000° higher. should pointed out however that 
Williams (12), the basis color temperatures, places the temperature 
the stars 15,000° while the Greenwich observers (6) place them even 
higher that means certain the present determination error. 

For the type stars may recalled that Petrie has derived temperatures 
28,000° 36,000° while Payne (10) earlier work, gives 35,000° for 
the stars class. The present determination (32,000° 40,000°) 


| 
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reasonably good agreement with these values and, since the discovery the 
binary character certain Wolf Rayet stars has explained the high values for 
the (65,000°) given Kuiper, now becomes possible speak with some 
confidence the order magnitude the temperatures this part the 


stellar sequence. 
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RADIO METHOD STUDYING THE YAW SHELLS' 


Abstract 


new method investigating the angular motion shell flight 
described. This radio method wherein the fuse the shell replaced 
high frequency transmitter and the radiated signal recorded ground 
station. The theory the method given, the equipment described, and the 
procedure the analysis the records discussed. Results experiments 
are quoted and values the force coefficients are deduced. addition those 
coefficients usually determined jump card experiments, certain other coeffi- 
cients, involving spin, are evaluated. The radio method measuring yaw 
compared with other methods. 


Introduction 


This paper contains account certain work ballistics carried out 
the authors the Department Physics, University Toronto. 


The subject the work new method investigating the angular motion 
shell flight. This method was made possible techniques originating 
the radio proximity fuse development. Early the fuse development 
was found that the spin the shell could determined, and considerable 
attention was given this both the U.S. and work was 
done the University Toronto effort correlate certain features 
spin recordings with secondary oscillations the shell. The essential feature 
the present method, however, was put forward Mr. Thomson 
the Directorate Naval Ordnance (U.K.), who suggested method for the 
measurement the yaw angle shell. 

the radio method investigating shell’s flight the fuse the shell 
replaced high frequency radio transmitter and the radiated signal 
recorded suitably located ground station. The angular oscillations 
the shell produce corresponding changes signal level, and this paper describes 
how these fluctuations signal strength may used determine the angular 
motion. theory the radio method given, and, after certain sim- 
plifying assumptions are made, formula suitable for application experi- 
mental results derived. 

The angular motion shell flight has been the subject considerable 
study the past. Most previous work has, course, been devoted prob- 
lems immediate practical concern the design projectiles. 
fundamental study was made Fowler, Gallop, Lock, and Richmond (2) 
who published their results just after the Great War. their paper, system 
forces was assigned, the equations motion were set and solved approxi- 
mately, and the theories were tested experiment. this paper the same 
‘force system and angular motion assumed, and there follows summary 
those results that are employed the radio method. 


received original form September 1947, and, revised, February 11, 1948. 
Present address: Institute Advanced Studies, Princeton, 
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II. The Theoretical Background 
THE SYSTEM 


Let the instantaneous direction the axis the shell represented 
Fig. OA, and the instantaneous direction motion the center 


(i) The plane AOP the plane yaw; 
(ii) The angle AOP the angle yaw, denoted the letter 
(iii) The angle rotation the plane yaw about the azimuth, 
denoted the letter 


Fic. 


For the moment, regard the shell rest and the air streaming past 
with uniform velocity, equal and opposite that the shell. The aero- 
dynamic forces can represented through symmetry single force 
located the plane yaw, and acting through the center pressure, 
the axis OA. The force can replaced, shown Fig. 

(i) The drag which the direction the air stream; 
(ii) The cross-wind force which acts through the plane yaw, 
perpendicular and positive when tends move the 
(iii) The moment about perpendicular the plane yaw, and positive 
when tends increase 


The following forms are assumed for and 


where the density the air, the radius the shell, and the velocity 
sound the undisturbed medium. 


These equations are the most natural form make the force coefficients 
through symmetry, the factor sin included explicitly order that 
may have non-zero limits discussion the dependence 
the force coefficients the constants the medium (diameter molecules, 
ratio specific heats, etc.) given 1.01 the paper Fowler al. 


Pp 
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practice the direction the axis the shell changes relative the 
direction motion. that and are unaltered the angular 
velocity the axis, but that the effect the angular motion the axis the 
shell can represented the introduction damping couple, called 
the moment due and acting such way directly 
diminish the resultant angular velocity, the axis the shell, shown 


Fig. 


There possible effect due the spin the shell about its 
assumed that the spin does not affect any the preceding force components 
but gives rise additional components. There will couple which 
tends reduce When the shell there will sideways force K’, 
analogous that producing swerve tennis ball, which need not act 
through but which must act normal the plane yaw, not alter 
through and couple about The complete effects spin can repre- 
sented, therefore, the addition the couples and and the force 
acting shown Fig. 

These forces are represented the equations: 


(2) 
before, the force coefficients are physical dimensions. Moreover, these 
force coefficients may depend appreciably other arguments, such 
and 
deriving the equations motion, Fowler al. were able simplify 
the algebra defining new coefficients the equations: 


ANysin 


where and are, respectively, the moments inertia the shell about 
longitudal and transverse axes through the center gravity, and the 
mass the shell. 
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The complete force system seen consist seven components; these, 
far the most important are and then, somewhat less important, 
Only one relation between the forces available, that between and 
For shells given shape moving given manner, the forces and are 
independent the position, the center gravity, while the moment 


where positive when nearer the base than For small 
this equation reduces to: 


This the only method finding 


fr) 


DESCRIPTION THE ANGULAR MOTION 


consider here only what Fowler al. call Type 
type motion occurs the early part all trajectories, and for the whole 
trajectory whose initial elevation less than 45° and muzzle velocity 
greater than few hundred feet per second. The characteristic feature 
Type motion that the axis the shell and the direction motion the 
deviate only small angles from the tangent the corresponding plane 
trajectory. simplify the treatment employed here, the assumption 


made that the directions motion the plane trajectory and the actual 
trajectory are parallel for that region the trajectory under investigation 
(up six seconds from the instant firing). This assumption can justified 
experimental accuracy this case reference Sections 3.64 and 3.65 
Reference The following system axes chosen; the axis the 
direction motion the center gravity; the axis the upward normal; 
and the axis horizontal and the right viewed from the gun (Fig. 3). 


Fic. 


the axis the shell, which represented the unit vector has 
components and (the direction cosines OA). The angle yaw is, 
definition, the angle between and OA; and the angle orientation, 
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specified the angle between the plane and the plane 102. and 
are related and by: 


The assumptions made Fowler, Gallop, Lock, and Richmond treating 
motion Type are the following: 
(i) The axis makes slowly varying angle with the horizontal, 
small compared the spin 
(ii) The yaw small (less than, say, 7°), 
(iii) The total effect the components I’, small, that all terms 
other than those the lowest order can neglected. 
The stability factor for spinning shell defined 
(7) 
found that the condition for stability that greater than unity. 
This is, course, analogous the condition for stability spinning top. 
When this condition satisfied, the following solution found for the 
equations motion the terms the direction cosines, and (these 
equations follow directly from those 4.0 the Fowler paper) 


t t 


where: 
0 


(9) 


Bo, and are arbitrary constants integration, and 


The subscript zero used give variable its value muzzle. 
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Equation (8) seen that the vertical component and the horizontal 
component the motion the axis, are composed two periodic terms. 
One, called the nutational term, consists damped oscillation instan- 
taneous frequency and instantaneous amplitude The 
other, called the precessional term, damped oscillation frequency 
and amplitude B). The ratio these two frequencies for 
the projectiles employed here was about exaggerated graph 
shown Fig. with the contribution each term shown separately. 


term. term. 


Fic. 


deduce all the force coefficients previously defined except fr, fr, and 
can derived from m2) and m2); may determined from 
the values found for various positions the centers gravity the 
shells. and may determined evaluating the exponents decay 
chronograph measurements the retardation the forward motion the 


shells. 


III. The Radio Method 
QUALITATIVE DESCRIPTION THE METHOD 


The Thomson radio method yaw measurement depends the fact that 
for radio fused shell, the radiation direction making angle with 
the axis the shell very closely proportional sin This most nearly 
true the region 10° from the axis the rear. 

The electric axis the shell radiator comes close to, but coin- 
cides with, the physical axis symmetry. small angle between the two 
will assumed deriving the fundamental equations, but will taken 


zero the present section. 

receiver located the instantaneous axis the shell extended, 
signal will received. For small angular displacement the axis, the 
electric field intensity the receiver will proportional the sine the 
angle between the axis and the line joining the shell and the receiver (line 
sight). the receiving antenna oriented respond, say, the vertical 
component the electric field, the signal received will proportional the 
vertical component the angular displacement (Fig. Thus the received 
signal will zero when the axis deflected horizontally, the right left, 
and will reach maximum value when the axis deflected vertically either 


172 
m™m. 
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There method distinguishing between these and down 
deflections, and there means determining the absolute value the 
angular displacement from receiver record. 


Line Sight 
Vertical 
Component 
Receiving 
Polar Pattern 
Fic. 


These difficulties are overcome for yaw measurement locating the 
antenna that there always angle between the axis and the line 
sight excess the expected yaw angle, and orienting the antenna 
respond only radiation polarized the plane this angle. this, 
the antenna located the plane the trajectory, distance ahead the 
gun, and oriented respond only vertically polarized radiation, Fig. 


Receiving Antenna 


The signal will then proportional m), where the angle 
between the line sight and the tangent the trajectory, and the vertical 
component the angle between the axis the shell and the tangent the 
trajectory. The mean signal level approximately 
Hence the ratio equal the ratio between the instan- 
taneous and mean signal levels (Fig. From knowledge the trajectory 
and the geometry the receiver and gun locations, possible find the 
angle the shell moves through space. Then knowing the value 
possible determine proper choice Q.E. and the distance gun 
antenna, possible arrange that stays reasonably constant over 


Axis 
2 
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desired region the early part the trajectory. For example, certain 
the experiments discussed later, lay between and 3°, from one five 


seconds time flight. 


Fic. Receiver record. 
DERIVATION THE FUNDAMENTAL EQUATIONS THE RADIO METHOD 


(a) Assumptions 

The following are the assumptions the derivation the fundamental 
equations the radio method: 

(i) The polar pattern the radiation from the shell sinusoidal, i.e., the 
field strength proportional the sine the angle between the electric axis 
the shell radiator, and the direction propagation. 

(ii) The electric axis the shell inclined the geometric axis, 

(iii) The radiation from the shell plane polarized; i.e., lies the plane 
containing the electric axis and the direction propagation. 


(iv) The receiving antenna responds only the component lying 
the plane containing the antenna, and the direction propagation. 


(b) Antenna Voltage Due toa Shell Space 

For the general case radiating shell, located and oriented arbitrarily 
space, define the following vectors: 

(i) unit vector coinciding with the geometric axis the shell, 

(ii) unit vector coinciding with the electric axis the shell, 

(iii) unit vector the direction propagation, i.e., line joining the 
C.G. the shell and the reference point the receiving antenna, 

(iv) the electric field vector the receiving antenna, given by: 


(10) 


where scalar proportionality factor depending power, wave length, 
etc., and the distance from the shell the receiving antenna. 

(v) unit vector the direction the antenna axis. The voltage 
induced field infinitesimal dipole can expressed as: 


+ 

. 
Signal Zero 
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For 


where scalar quantity depending wave length, units, etc. 
directional array, may write: 


where the polar angle measured from 

the azimuthal angle measured from arbitrary plane through 
scalar function, such that represents the gain 
the array relative the dipole. 

Then, substituting for from Equation (10): 


(11) 


determine the components these vectors employ the system 
axes, 01, and previously defined (Fig. 8). Let O1, represent 


unit vectors their axes. 


makes angle (the yaw angle) with O1, and the plane and 
makes angle (the angle orientation) with the 102 plane. 


the plane and perpendicular and Ob. The electric 
axis passes through the center gravity, and makes angle with 
Owing the spin the shell, rotates with angular velocity about 

Then the triad, reckoning the angle from the plane: 

transform the vector into the 123 triad require the components 

Oa, Ob, Oc, referred the 123 triad: 


‘ 
\ = \ a 
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and Ob, are determined unit vectors perpendicular small, 


Substituting the above values for Oa, Ob, and Oc, Equation (13), with 

lies the 102 plane, and the angle between and and therefore: 


lies the 102 plane, inclined the angle that 


Substituting these values Equation (11), have: 


arrange the experiments that remains small, that sin 
cos and compressing sin into (since nearly constant 
for actual trajectory), this may written: 


When the spin not immediate interest put equal zero obtain: 


THE EQUIPMENT 


The equipment required for measurement yaw the radio method 
the following: 

(i) radio fuse that can fitted standard artillery shells and that 

(ii) receiving antenna, responding vertically polarized radiation, 

(iii) radio receiver, tuned the frequency the signal from the fuse, 
and whose output voltage proportional the input voltage from the antenna, 

(iv) cathode ray oscilloscope, permitting the output voltage photo- 
graphed, 

(v) timing device, indicating the film record (1) the time when the 
gun fired and (2) the time any subsequent point the film. 

The equipment developed for this work described this 
Part (a) the fuse, radio sonde, discussed, the antenna, and (c) the 
receiver, and its associated apparatus. 


Fic. yaw fuse. 


PLATE 


‘J 


7 
Z 


i 
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(a) The Radio Sonde 

The requirements the radio sonde are: 

(1) The polar radiation pattern the rear shell must sinusoidal, and 
must figure revolution whose axis nearly coincides with the axis 
the shell. 

(2) Power radiated must sufficient give signals well above the noise 
level the receiver, even when the shell the largest distance from 
the receiver. 

(3) The frequency must capable adjustment the laboratory, that, 
firing, the sonde frequency will lie nearly possible the mid- 
band the receiver. 

(4) The frequency drift the unit over the useful time flight must 
small. 

(5) The variation output power over the useful time flight must 
small. 

photograph the sonde shown Fig. oscillator circuit 
conventional, but course employs rugged (V.T. Fuse) components. The 
battery regular V.T. Fuse battery. 

Requirement (1) was attained using the lowest possible frequency. 
The portion the polar pattern actually used that not exceeding 20° from 
the axis the shell the tests were made large open 
field check this region. The shell was supported polystyrene rod about 
The axis the shell was horizontal and the shell was rotated about 
vertical axis. Relative field strengths were then measured few hundred 
feet away. tests indicated that the polar pattern was sinusoidal 
the required region, closely addition, within the accuracy 
these tests, the electric axis and the geometric axis the shell appeared 
coincide. discussed Section 5(a), this not strictly accurate. 

Requirement (3) was met taking special care with the oscillator coil, 
and using, for final tuning, number nose caps different sizes. This 
gave sufficient adjustment bring all the frequencies within the desired range. 

Requirement (4) small frequency drift was checked following number 
units flight with narrow-band receiver. tests did not reveal 
any frequency drift. 

Requirement (5) power drift was checked the laboratory. Over 
relatively long periods (several minutes) change power output could 
detected. 
(b) The Receiving Antenna 

The requirements the receiving antenna are: 

(1) directivity pattern that discriminates azimuth against reflections 
from behind and from the sides. 


(2) Independence effects due variations ground conditions. 
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(3) the region maximum reception response that independent 
elevation over range about 10° and azimuth over range 
least 4°. 

(4) radiation resistance match the ohm characteristic impedance 
the coaxial cable the receiver. 


(5) Maximum possible gain over half-wave dipole. 
(6) Portability and ease assembly the field. 


(7) Provision for tilting the vertical plane investigate any part the 

trajectory. 

was decided use 90° corner array, but provide with ground 
screen order ensure the second requirement above. This type antenna 
was first investigated Kraus (4), whose conclusions can summarized 
follows: (1) the physical dimensions the array are not critical, (2) power 
gain possible using 90° corner, and (3) for dimensions smaller than 
two wave lengths, the performance corner arrays equals that parabolic 
reflectors and horns. 


The directivity patterns the antenna azimuth and elevation are 
shown Figs. and 11. 


Fic. 10. pattern antenna azimuth 90°). 


40° 


Fic. 11. pattern antenna elevation 0°). 


20° 
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(c) The Receiving 
The requirements the receiver are: 

(1) Ability find the signal within 0.5 sec. firing the shell. 

(2) sensitivity adequate for good photography for the trajectory region 
under investigation. 

(3) means recording the amplitude the radio-frequency signal. 

(4) Linearity output with input. 

(5) Insensitivity line voltage fluctuations. 

(6) method providing accurate time scale from the instant the 
gun fired. 

Fig. shows the block diagram the receiving system. 


Antenna 


D.C. 
Amplifier Amplifier 

Power 

Supply 


12. 


Block diagram receiving system. 


Since the yaw largest very early flight, essential that there 
little delay possible obtaining usable signal. Moreover necessary 
hold this signal once found, since tuning the receiver would affect its output 
and make yaw measurements impossible. Therefore was decided that 
manual tuning during shell flight was undesirable, and that fixed-frequency, 
broad-band receiver would meet the requirements best. Because the sondes 
may show certain frequency spread there was included tuning condenser 
the local oscillator. Before each round fired the receiver tuned the 
laboratory measured frequency for that sonde. Nevertheless bandwith 
few megacycles desirable, because sharply tuned receiver would respond 
frequency drift change signal strength which could not distin- 
guished from changes due variations orientation and distance the 
shell. Automatic frequency control was considered means eliminating 
this difficulty, but was found unnecessary. 


Monitor 
Amplifier 
Recording Camera 
Timing 
Unit 
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Preliminary calculations and experience with the proximity fuse indicated 
that ensure good signals over sec. time flight would necessary 
employ high sensitivity receiver. This problem high sensitivity and 
wide band was exactly that which had been overcome the construction 
radar receivers, and was decided that would simplest modify such 
receiver. The frequency was adjusted cover the desired range, and 
switch was introduced, thus permitting the gain altered six fixed 
steps. For linearity response, diode detector was employed. 


Recording the carrier level the radio-frequency signal requires direct 
coupling between the detector and final indicator, the cathode ray oscillo- 
graph. direct-coupled amplifier consisting single 6AC7 was employed 
for this purpose. 

Sometimes advantage listen the spin modulation the sonde. 
For this there A.F. amplifier and speaker. 


The d-c. amplifier followed two in. cathode ray oscilloscopes, one 
with green tube for monitoring, the other blue for photographing. The no- 
signal voltage the C.R.O. plates that the d-c. amplifier plate, 
175 Therefore, the high voltage must specially insulated make 
centering and focusing possible. horizontal sweep voltage may applied 
the monitor desired. Because the vertical motion the film the record- 
ing tube rotated 90° that horizontal deflections correspond vertical 
deflections the monitor. 

order relate the instantaneous position the shell space the 
corresponding point the trace necessary put time scale the 
records. This was provided series timing marks intervals 1/100 
sec., about 1/6 in. apart the film. These are obtained photographing 
narrow slit illuminated neon bulb, pulsed from independent timing unit. 


The master timekeeper 1000 c.p.s. tuning fork oscillator (General 
Radio Type 813-A). The output the oscillator fed into peak clipper 
and thence step circuit. careful adjustment the cathode voltage 
the step tube, the step circuit made trigger asymmetrical multi- 
vibrator the beginning the seventh step. Another four milliseconds are 
allowed for the step circuit discharge; thus the full cycle completed 
1/100 sec. The multi-vibrator delivers 300 pulse duration 
the neon bulb. Examination the wave form the voltage across the 
neon bulb shows that the firing time repeats itself better than 


fiducial mark for the time scale obtained momentarily switching 
resistor parallel with the plate load the d-c. amplifier. This introduces 
sharp rise short duration the trace. The switch operated the 
pressure from the blast the gun. The moving part copper disk about 
in. diameter. The whole switch mounted box and placed the 
side and slightly behind the muzzle the gun. The time delay switch 
operation was measured separate chronograph trial and the necessary 
correction was made for each round. 


| 
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The zero signal position the C.R.O. spot extremely sensitive voltage 
fluctuations, especially the I.F. screen voltage, the d-c. amplifier heater 
and the high voltage the C.R.O. The primary source power for most 
the trials was cycle, 1500 w., Homelite gasoline generator. This has 
quite stable operation for its type, but surges occurred frequently. 
Therefore was decided use power supplies stabilized negative feedback. 
similar type unit used the 700 supply the and here 
also the stability more than adequate. The heater voltage the d-c. 
amplifier supplied storage battery. The heater voltages the 
R.F. and I.F. amplifiers are not critical. The over-all system such that rapid 
generator fluctuations not affect the spot position appreciably. slow 
change causes steady drift but this negligible over the sec. time flight 
single round. 

The Recording Camera employed the General Radio Oscillograph 
recorder, Type 651-AE, with lens 651-PI. This continuous drive camera 
taking either mm. film; mm. film (Kodak Super XX) was used 
throughout this work. For photography the lens opened f/1.5. 
The maximum writing speed about in. per sec. 


The complete equipment held shockproof steel rack and can 
mounted three ton Dodge truck which complete equipment for trials 


carried. 


IV. The Experiments and Experimental Results 


EXPERIMENTS 


DESCRIPTION THE 


Experiments were performed the radio method measure yaw 
the pr. shell. The firings were performed Nov. 15, Nov. 23, and Dec. 
1945. The Artillery Proof Establishment the Inspection Board the 
United Kingdom and Canada, Hamilton, Ontario, supplied guns, propellant, 
and gun crews. The firings were held the Winona Rifle Range the 
Department National Defence. 

The gun employed was the pr., firing H.E. MK. streamline shell. 
All rounds were fired with normal full charge quadrant elevation 20°, 
from barrel the first quarter wear. The gun was located about 300 
yards from the shore, and pointed over Lake Ontario. The antenna was 
located 227.4 ft. front of, and 6.1 ft. below, the muzzle the gun, the 
vertical plane containing the trajectory. mound earth, ft. high 
ft. thick was located between the antenna and recording station, and the gun. 


The velocities and retardations all rounds fired were measured officers 
the Canadian Armament Research and Development Establishment, 
Valcartier, Quebec, using the counter chronographs and photoelectric screens 

normally used these measurements. the trial, measurements were 
made temperature, relative humidity, and wind velocity ground level. 
The mass, the transverse and longitudinal moments inertia through the 
center gravity and B), and the position the center gravity, were 
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determined for all rounds order determine Equation (5) 
shells three different mass distributions were fired. one group sulphur 
filling was employed, another, lead, and the third, steel plug was 
placed the nose beneath the fuse. This produced range positions 
the centers gravity about in. 


DESCRIPTION THE TRACE 


typical trace reproduced Fig. following are salient points: 

(i) the region one three seconds, three frequency components are 
evident: 

(a) high frequency, about 250 c.p.s., with fairly small amplitude. 
This the spin 

medium frequency, about c.p.s., with amplitude the order 
one-quarter film This the term the vertical com- 
ponent, the shell’s angular motion (see Section 2). 

(c) low frequency, about c.p.s., with amplitude the same order 
(b) above. This the term. 

(ii) The component gradually decreasing amplitude, until 
three seconds, when has almost wholly disappeared. There is, meanwhile, 

(iii) the region one two seconds, the level the entire trace shifting 
downwards. Later, remains almost constant, when the decrease due 
nearly balances the increase due increasing 

Certain other points should noted. the middle the film there 
series small rectangles. These are the timing light marks, 1/100th sec. 
apart. They are numbered every 1/10th sec. Before the point marked zero, 
the trace the lower edge the film, steady straight line. zero 
time, and following for 0.15 sec. there are several jogs the film. These 
are caused the pressure switch, actuated the gun blast. about 0.2 
sec. the first signal received, but great send the trace off scale. 
about 0.8 sec. the signal has decreased and scale after that. The 
sudden jump the height the trace 2.2 sec. caused scale change. 
When the receiver operator saw the signal becoming too low the monitor 
oscilloscope, switched the receiver higher gain position, causing this 
discontinuity. 


PROCEDURE AND THEORY ANALYSIS 


The voltage the receiving antenna, due the shell space, is: 


The trajectory the shells, course, must known permit analyses 
the records. The trajectory data specifically required are: the distance 
along the line sight from the antenna the shell; the angle between the 


PLATE 


Fic. Photographic trace Round 67. 5.3 sec. 


uw 1.0 0.9 08 0.7 0.6 05 
4.\ 4.0 39 3.8 3.7 3.6 3.5 


III 


13, continued, 5.3 10.8 sec. 


1 
1 
6.4 82 81 692 1.9 7.8 7.7 
9.0 89 8.8 8.7 8.6 65 84 
1 
9.6 9.5 94 9.3 9.2 9.0 
02 10.0 2.9 2.6 9.7 
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line sight and the tangent the trajectory, and the angle between the 
line sight and the reference axis the antenna, (Fig. can 
determined from the trajectory and from the position the antenna relative 
the gun. From these data, and from the calculated response pattern 
the antenna, values and were tabulated and plotted against time. 
1 
and are the values and one second. was necessary cal- 
culate the trajectory these shells, since the shape the yaw fuse different 
from the fuses regularly employed. short-arcs method was used for this 
purpose, the ballistic coefficient being determined separate retardation 


trial. 
(a) Preliminary Steps 

facilitate measurements, enlarged copies the trace were made 
projecting its image continuous graph paper, and tracing pencil. 
The complete motions the oscilloscope spot and the timing dots were care- 
fully reproduced. magnification about seven was employed, and, since 
the distance between sprocket holes the film in., the working traces 
were in. wide. 

operation the graph paper was smooth out the spin 
tion, either estimating the mid-point directly drawing the envelopes 
the peaks and taking the curve midway between. 


(b) The Analysis 
the height the curve above the zero signal level, linear* function 
the input voltage 


Equation (19) 
m), spin being smoothed out. 


Hence 


useful define two new quantities, follows: 


(i) The variable where 


The factor which constant, will included hereafter, that 
1 


the dimensions become inches displacement the graph, per degree 
angular displacement the shell’s axis. That is, write 


There slight the input vs. output the receiver, but since this can eliminated 
addition small constant use the curve itself, the above discussion 


not error. 
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was calculated from for several hundred points each trace, using 


Equation (21), and the curve 


(ii) The variable where 
where constant, the first approximation 
The purpose the remaining preparatory work was replace the curve 


since the contribution very much smaller than the contribu- 
tion Km. isa more suitable graph for analysis the shell’s motion. 


was not difficult, inspecting the curve locate approximately 
the times the maxima and minima the low frequency modulation. 
was assumed, for this approximation, that the time midway between 
maximum and minimum, this low frequency gives contribution the 
trace amplitude. inspection, drawing the envelope the peaks 
the medium frequency was possible estimate the contri- 
bution this frequency component the time Thus good estimate 
could made the trace amplitude for zero vertical component yaw, 
But then Ky, and since known all times, this 
within the finally determined. 

From the values obtained series points the trace, from the 
tabulated values and from the values was calculated Equation 
versus time was then superimposed the original tracing 
Two curves, and the envelopes the maxima and minima, 
respectively, the term, were drawn. graph typifying the 
curves and the envelopes and shown Fig. 14. 


Ge 
Ge 
Time 


Fic. Gi, and vs. time. 


The equations these curves are: 
t 


a 
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(25) 


adding Equations (24) and (25) get: 


tz 


(c) The Low Frequency 


tmax 


Combining Equations (27) and (28), 


tmin 


fmax 


where the time maximum and the time the adjacent 


minimum. Assuming m2) constant over this interval, may write: 
2(tmin 


This can, course, applied equally well from minimum toa maximum 


Period 
Low Freq. Secs. 


Time, Secs. 
Fic. 15. low frequency vs. time. 


Round 89. 


The values for typical round (No. 89) are plotted against 
time Fig. 15, and smooth curve drawn through them. The value 
employed determine force coefficients were taken from such 


-250 2S 
° 
3.0 4.0 


186 CANADIAN JOURNAL RESEARCH. VOL. 26, SEC. 


curves for each round. These values, and their corresponding times, are 
shown Table for each round. 


TABLE 
oun oO. N2), C.p.S. ny Nz), C.p.s. m2), 

22.01 2.31 1.60 

22.01 2.46 1.18 

22.42 2.19 1.60 

22.35 2.24 1.60 
20.41 2.89 0.98 
20.37 2.50 1.35 
20.75 2.50 
20.55 2.42 1.60 
20.73 2.25 1.60 
21.19 2.19 1.60 


(d) The High Frequency 


was observed, for all rounds analyzed, that m2) was constant over 
the region which could observed the trace. Thus, was readily 
determined noting the times its minima and maxima, and thereby cal- 
culating the mean period from the time taken for complete cycles. 
The range generally used was that from one point low frequency cycle 
the corresponding point the next cycle. 


Table the values are shown for each round. 


(e) Determination 
From knowledge 


t 
ar | (m = dt 


ty 


can determine the times which 


Then from Equation (26) such times: 


(30) 


This determines differs from more than 5%, the variation 
with time may have affected sufficiently impair the accuracy the 
Equation (30) must regarded second approximation the true 
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The process must then repeated. almost every case, repetition was 
found unnecessary. The values prior the first scale change are shown 
Table II, and with them, the corresponding K’. 


TABLE 


inches per degree inches per degree 


Round No. 


0.90 


(f) Determination the Amplitudes 
Adding and subtracting Equations (24) and (25) have: 


Equations (31) and (32) were used determine the yaw amplitudes. 


semilogarithmic scales Fig. corresponding least squares straight 


te 
round are plotted Fig. 16. 


THE DERIVED QUANTITIES 


The quantities determined direct measurement from the trace are the 


addition, the following quantities were determined, either laboratory 
measurement, independent measurement the trial: 
The moments inertia, and the shell, 
The position the center gravity, 
The velocity and mean ballistic coefficient the shell, 
Meteorological data the time the trial. 


0.92 

1.22 

2.18 2.24 

0.96 0.91 

0.86 0.85 

0.72 0.70 

1.47 

1.04 
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Certain these measured values are shown Table 


TABLE III 
Position center 
mass Axial Transverse 
R a Mass moment moment one 
t 
Normal 6.64 24.97 40.42 418.9 1487 
6.64 24.98 40.19 418.7 1486 
6.62 24.95 39.93 413.0 1499 
6.62 24.94 40.32 415.6 1491 
Back 6.31 25.00 40.10 441.7 1495 
6.33 24.97 40.21 445.2 1486 
6.25 24.94 437.9 1498 
6.27 24.98 39.94 439.1 1488 
Forward 7.03 25.08 39.99 440.1 1490 
7.03 25.11 39.94 433.6 1490 


The velocities the shells, tabulated functions time, were calculated 
the standard methods external ballistics. The force factors were then 
determined from the above data, described Sections and 

(a) The Overturning Moment Coefficient 


From Equation (9) 


and since 
fi 


may obtained independently the relation using the 


measured values and the spin Using determined this way, 
and the measured values m2), may evaluated, since 
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evaluation and The values and found both methods 
are given Table IV. 


TABLE 


Round From From From 


1.60 24.32 24.12 13.61 12.41 6.64 
1.18 24.47 24.31 13.85 12.93 6.64 
1.60 24.61 24.61 12.90 12.94 6.62 
1.60 24.59 24.60 13.36 13.41 6.62 
1.35 22.87 22.96 14.55 14.55 6.33 
23.33 13.93 14.33 6.25 
1.60 22.97 22.87 14.01 13.40 6.27 
1.60 22.98 13.20 13.00 7.03 
1.60 23.38 23.28 12.94 12.33 7.03 


(b) The Cross-wind Force Coefficient 

comparing the values for shells three different positions the 
centers gravity, the values may deduced from Equation 
the data Table was plotted against position center gravity and 
straight line fitted least squares (Fig. the slope this straight 
line was found that 2.0, whence (as defined equation (3) 0.18. 


Moment Coefficient fm. 


Inch 
Dist. Mass Centre from Base 


Fic. vs. distance mass center from base. 


(c) The Damping Factors 
The quantities are called the ‘damping factors’, and are related 
their corresponding force factors Equation (3). 
The determination has been described above. determined from 


14 
12 
6.0 6.2 6.4 6.6 6.8 70 72 
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The factors and may determined from the observed values 
and using Equations (9). The best exponential curves the 


form were fitted the observed values least squares. This 


involves neglecting the effect variations and assuming that 


the exponents are linear functions time. Consideration the rate 


change shows that changes only fraction per cent 


per second. The assumption regarding the exponents can justified 
regarding the values obtained averages. The semilogarithmic graph 
against time Fig. (17) shows that this procedure reasonable. 

hand 

exponents, are tabulated. The values and are also shown. 

will noted that negative. The significance this that the 
force acts behind the center gravity, rather than front. 


SPIN 
Modulation Amplitude 


Equation (18) contains term cos which component varying 
the spin frequency. For those times when cos follows that: 
where the height the spin modulation the trace and height 
trace when averaged over spin cycle. Measurements the traces show 
that, within experimental accuracy, remains constant throughout flight, 
that the spin term the signal adequately explained postulating the 


existence angle between the electric axis and the spin axis the 
shell. For Round 67, 0.27°, and for Round was 0.19°. 


(b) Measurement Spin and Spin Decay 

The small amplitude spin modulation and the relatively low film speed 
make statistical method measurement necessary. Every tenth spin 
cycle the trace was numbered, and the observed time for cycles, 
was obtained. was thus calculated, and Fig. shows typical graph 
spin against time. For each round curve the form was fitted 
determining the least squares straight line: 


(34) 


The values and found are listed Table most rounds 
this equation was applied over five six second interval. 
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TABLE 
Round No, rev./sec. Round rev./sec. 
0.0092 257.0 0.0109 
256.2 0.0095 257.2 0.0109 
259.0 0.0107 0.0123 
254.8 0.0108 259.3 0.0115 
257.8 0.0103 254.9 0.0087 
250 
245 
240 
a 
> 
c. 
230 
225 
220 


Fic. Round 12. 


Application Equation (35) over shorter time intervals showed that there 
was significant trend This trend was investigated determining 
for two-second intervals, five rounds. About points were 
all but few the range sec., and smooth curve drawn through 
them. The mean distance the points from the curve was about 3%. 
Using the relation connecting and the spin reducing couple, 


the curve versus time Fig. was obtained. 


(c) Skin Friction 


Spin decay results from friction due the viscosity the air and 
possible estimate from the theory skin friction. The usual hydro- 
dynamic treatment terms the boundary layer theory (3) according 
which the shell surrounded thin layer air which the inner surface 
rest relative the shell, and the outer surface moves with the velocity 


W 


21s 
Time, sec.— 
Vi 
al 
al 
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TABLE VII 


CALCULATION 


Time, 
sec. 


Mach number V/a, velocity sound the undisturbed medium and 
velocity air stream. 


(36) 


the air stream. Although the viscosity small, the normal velocity 
gradient may large, and the integrated effect the shear stresses can give 
rise appreciable force. The drag written the form: 


(37) D= ’ 
where dimensionless coefficient, 
effective length shell length fuse correction). 
From this follows that the spin reducing couple 


and comparing this with Equation (2): 


(38) 


For Falkner has modified theory Karman’s show that 
(39) 


where the Reynold’s number, 


The values calculated from Equation (38) and those from the measured 
values are shown Table VII and Fig. seen that the magnitudes 
are reasonable agreement but the shapes the curves are 
should roted that the above theory for the drag does not take into account 
any possible changes the nature the equations flow that occur the 


Time, V/a from from from from 
sec. drag drag 
0766 0.898 0.0715 0.0837 
0779 0.882 0.0719 0.0839 
0787 0.866 0.0733 0.0842 
0797 0.852 0.0754 0.0842 
0809 0.841 0.0772 0.0844 
0814 0.829 0.0801 0.0847 
0819 0.820 0.0837 0.0850 
0824 0.810 0.0872 0.0850 
0832 0.802 0.0914 0.0850 
0834 
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trans-sonic region. The minimum 6.5 sec. corresponds Mach number 
about one. 


Fic. 20. Coefficient spin couple. 


interest see how much the drag due skin resistance contributes 
the total drag shell. The total force due mainly wave drag and 
given the empirical equation: 


where 
ballistic coefficient, 


and 


velocity sound under standard conditions). 


From Equation (2) follows that; 
Totaldrag 


found that the skin drag contributes from 25% the total drag, 
the fraction being least 


(41) 


Summary and Conclusions 
SUMMARY PRECEDING RESULTS 


The system forces acting shell flight described, after the analysis 
Fowler, Gallop, Lock, and Richmond. The solutions the equations 
motion found those authors are quoted, indicating the manner which the 
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forces affect the angular motion. radio method measuring this motion 
over wide range distances from the gun described, the theory the 
method derived, and the equipment developed for the measurements 
described. account certain experiments using this method given, 
and the photographic record described. The theory and examples the 
analysis are given. Results for some rounds are quoted. From these 
results, the force factors and damping coefficient have been deduced. The 
overturning moment coefficient has been determined velocity 1350 
ft. per sec.; and from values for shells various mass distributions, the 
cross-wind force coefficient has been found. Average values the 
damping factors, and have been deduced evaluating the exponents 
the decay the yaw angle. accurate measurements spin, the spin 
decay exponent has been determined, and the spin decay coefficient, has 
been evaluated different times. has also been calculated from hydro- 
dynamic theory, and comparison made with the values determined 


spin decay. 


COMPARISON THE RADIO AND OTHER METHODS 


Prior this development, two other methods yaw measurements have 
been employed. These are the jump card method and the photographic 
method. The jump card method consists firing shell through series 
cardboards, spaced down range. measurements the shape and 
orientation the resulting holes the cards, possible determine 
and The photographic method consists firing shell down range, 
between pairs cameras which photograph the shell perpendicular planes, 
number points down the range. From these photographs, and 
may determined. 

possible this point compare the radio method with the jump 
card and photographic methods yaw study. 

The relative advantages the radio method are: 

(a) permits measurement over the whole trajectory; 

provides continuous record; 

(c) mechanical contact with the shell required; 

(d) permits measurement the spin and spin decay; 

(e) The accuracy all measurements comparable with either card 
photographic method, and for number quantities, e.g., nutation 
and precession frequencies, damping factors, better; 

(f) Measurements are not confined flat trajectory. 

The relative disadvantages the radio method are: 
(a) requires special fuses and equipment; 
(6) The labor and time required for analysis considerably greater; 
(c) Analysis requires independent determination trajectory data. 
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Abstract 


Free air, Bouguer, and isostatic anomalies are listed table principal 
facts for stations established the Dominion Observatory the area 
between latitudes and degrees, bounded the east Hudson Bay and 
the west Athabaska, Slave, and Mackenzie rivers. There some evidence 
that the predominantly negative anomalies stations the Pre-Cambrian may 
due lack isostatic adjustment following removal glacial loads. 


Introduction 


The first gravity observations northern Canada were made the 
Dominion Observatory 1921 when series pendulum stations were 
established along the water routes between Edmonton and the Mackenzie 
River delta Observations were made seven stations the Arctic 
far north Cambridge Bay Victoria Island combined gravity and 
magnetic airborne expedition 1945. Three determinations were made 
points along the Hudson Bay railroad 1946. Altogether stations have 
been established the area between latitudes and degrees, bounded 
the east Hudson Bay and the west Athabaska, Slave, and Mackenzie 
rivers. The results all observations are listed Table and are shown 
Fig. 

interpret the results gravity observations the value gravity 
compared with theoretical values computed for the latitude and elevation 
the station. Several methods reduction may employed. observed 
value minus the computed value the gravity positive anomaly 
represents excess gravity and negative anomaly deficiency. this 
report three methods reduction have been used, the free air, Bouguer, and 
the Hayford and Bowie isostatic reduction for depth compensation 
113.7 km. theoretical value gravity sea level for each station 
computed from the International Gravity Formula* and shown the 
sixth column Table The corrections that are applied this value 
obtain the computed value gravity the the station are 
shown the next five columns the table. free air computed value 
obtained simply applying negative correction for the elevation the 
station above sea level account for the diminution gravity due the 


received December 12, 1947. 

Contribution from the Dominion Observatory, Ottawa, Canada. Vol. No. 
butions from the Dominion Observatory. Published permission the Director, Mines, Forests, 
and Services Branch, Department Mines and Resources. 


Theoretical gravity centimeters per second per second) sea level for each 
station obtained substitution the latitude the station the formula 978.049 
0.0052884 sin? 0.0000059 2¢), adopted the International Union Geodesy 
and Geophysics. this formula 978 .049 the most probable value gravity the equator, 
deduced numerous gravity observations many places over the earth. The coefficients 
sin® and are deduced from theoretical expressions for gravity the adopted spheroid. 
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increased distance from the center the earth. For the Bouguer (modified) 
reduction further positive correction made for the attraction the 
topography out distance 103.5 miles (Hayford Zone from the 
station. For the Hayford and Bowie (6) reduction the theoretical value 
sea level first corrected for the elevation the station and then for the 
topography the whole earth and finally for the compensation corresponding 
the topography. Bullard’s (1) method for computing the effects 
corresponding the lettered zones was used and maps Canada based 
Heiskanen’s charts for Zones were employed obtain the corrections 
for these more distant zones. 


Discussion Results 


perfect isostatic equilibrium prevailed and the vicinity the station 
there were large masses differing density from 2.67 (the value adopted 
the mean density the solid part the earth) the isostatic anomalies 
would tend negligible. Isostatic anomalies for the stations, grouped 
according the geological formations (on which they are situated), are shown 
Table are seven Palaeozoic stations, five with negative anomalies 


TABLE 


TABLE SHOWING ANOMALIES FOR THE VARIOUS STATIONS GROUPED ACCORDING 
TO THEIR GEOLOGICAL FORMATIONS 


Pre-Cambrian Palaeozoic Mesozoic 
Station number Station number Station number 
175 170 .042 +.001 
167 .004 +.005 .008 
Mean anomaly Mean anomaly .009 Mean anomaly 
with regard with regard with regard 
sign sign sign 
Mean anomaly Mean anomaly Mean anomaly 
without regard without regard without regard 
sign sign sign 


and two positive. The mean anomaly with regard sign —0.009 cm. 
per Five out seven Mesozoic stations have positive anomalies. Only 
one Mesozoic station has negative anomaly. The mean anomaly with 
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regard sign for the Mesozoic stations 0.002 cm. per sec.? All the 
nine Pre-Cambrian stations have negative anomalies and the mean anomaly 

Similar results have previously been observed Canada over all these 
formations, including the Pre-Cambrian. Pre-Cambrian rocks are the 
average most likely considerably heavier than normal, persistent negative 
anomalies over such extensive Pre-Cambrian area point departure 
from isostatic equilibrium overcompensation. 


The center the Pre-Cambrian area question may have been near the 
center the great Keewatin ice sheet during the Pleistocene period North 
America. Coleman (2) describes the Keewatin having different times 
three main centers, the first the west Chesterfield Inlet, the second near 
Dubawnt Lake, and the third near the mouth the Nelson River, the ice 
centers reaching estimated depth 12,000 has been well established 
geological studies that, owing the tremendous loads, glaciated areas were 
depressed far below present levels and that the land surface still stands lower 
than preglacial times Although the few statistics rate regional 
uplift that are available have been the subject debate (3) generally 
agreed that equilibrium will reached only after the termination thousands 
years. 

Lack isostatic equilibrium should then prevail regions postglacial 
upwarpings the crust, gravity anomalies would negative, and increas- 
ingly negative toward the area where the load was the greatest The 
gravity measurements made over this large area exhibit some extent this 
expected deficiency gravity. Churchill and Cambridge Bay, both stations 
near the former center the ice sheets, have negative anomalies —0.051 
cm. per and —0.042 cm. per respectively. the other hand 
Arctic Red River and Good Hope near the Mackenzie River delta have 
positive anomalies, which agreement with geological evidence that this 
region was the outer limits the ice cap and was not disturbed glacia- 
tion There also some evidence algebraic decrease the anomalies 
passing from west east. The area covered the Prairie Provinces 
the south and extending through the central United States also for the most 
part one positive anomaly. 

Conclusion 


The question whether the negative anomalies predominant 
Pre-Cambrian area are due lack isostatic adjustment following removal 
glacial loads cannot yet definitely answered. Although the existing 
measurements gravity suggest affirmative answer, greater density 
gravity stations, especially regions the northeast, may present different 
picture. hoped the near future make use gravimeter extend 
the system gravity measurements throughout this area such way 
answer some the questions raised this paper. 
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